AD-A116  383  TEXAS  UNXV  AT  AUSTIN  OCPT  OP  ELECTRICAL  ENGINEERING 
AN  ISOLATOR  IN  IMAGE  GUIDE. (U) 

MAY  82  N  CAMILLERIf  T  ITOH 
UNCLASSIFIED  UT-MW-N-a2-l 


N00014-79-C-0553 

NL 


Unc  Lassif led 


SF.CURI'^Y  CLASSIFICATION  OF  THIS  PAGE  (Whmn  Dmlm  Bnimnd) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1  REPORT  number 

LIT-MW-N-82-1 

2.  GOVT  ACCESSION  NO. 

F<D-.4|  l(p  Zi, 

3.  RECIPIENT’S  CATALOG  NUMBER 

5 

4  title  ("And  Subllllf) 

An  Isolator  in  Image  Guide 

5.  TYPE  OF  REPORT  &  PERIOD  COVERED 

Technical  Report 

8.  PERFORMING  ORG.  REPORT  NUMBER 

7.  author^*) 

N.  Camilleri  and  T.  Itoh 

8.  CONTRACT  OR  grant  NUMBERr«J 

N00014-79-C-0553 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

University  of  Texas  at  Austin 

P.O.  Box  7728 

Austin,  TX  78712 

10.  PROGRAM  ELEMENT,  PROJECT.  TASK 
AREA  S  WORK  UNIT  NUMBERS 

1  1.  controlling  OFFICE  NAME  AND  ADDRESS 

Office  of  Naval  Research 

Arlington,  VA  22217 

12.  REPORT  DATE 

11  May  82 

13.  NUMBER  OF  PAGES 

125 

14.  monitoring  agency  name  6  AOORESSCI/ dlMarcnl  ttom  Conitollini  Olticm) 

IS.  SECURITY  CLASS,  (el  thie  report) 

Unclassif ied 

I5a.  DECLASSIFICATION/ downgrading 

schedule 

16.  distribution  STATEMENT  (ol  (hl4  Report) 

Unlimited 

17.  distribution  statement  cof  (A*  •6«fr«c<  /n  Sloe*  20,  // d///»r«n/ from  RapwO  fc  !!  !;  f 

JUL2  19821 


18.  SUPPLEMENTARY  NOTES 


^‘1 

S  « 


19.  KEY  WORDS  (Cont/nu0  on  rovmroo  oido  If  nocooomrr  «ntf  fdonttfy  by  block  numbor) 

isolators,  image  guides,  Galerkin's  Method,  leaky  wave  filters 


50.  ABSTRACT  (Conllnum  on  rovoroo  otdo  If  noeoooofT  Idontlty  by  block  nunbor^ 

An  image  guide  consisting  of  a  layer  of  dielectric  and  a  layer  of  ferrite 
with  periodic  metal  strips  between  these  layers  is  analyzed.  Non-reciprocity 
of  the  phase  constants  caused  by  the  gyrotropic  nature  of  the  ferrite  makes 
it  possible  to  create  non-reciprocal  leakage  phenomena  by  means  of  the  metal 
strips.  These  phenomena  can  be  used  to  develop  a  distributed  isolator. 
Numerical  results  and  data  from  the  experimental  investigation  performed  are 
presented .  ^ 


DD  ,  1473  EDITION  OF  I  NOV  •(  IS  OMOLETE 

S/N  0102-014-6601  I 


UNCLASSIFIED 

SECUHITV  Ct  AMIFICATION  OF  THIS  PAOE  0»l« 


4 


MICROWAVE  LABORATORY  REPORT  NO.  82-1 


AN  ISOLATOR  IN  IMAGE  GUIDE 


by 

N.  Camllleri 
and 

T.  Itoh 


May  1982 


Sponsored  by  the  Office  of  Naval  Research 
Contract  No.  N00014-79-C-t)553 


Microwave  Laboratory 
Department  of  Electrical  Engineering 
Tlie  University  of  Texas  at  Austin 
Austin,  Texas  78712 


AN  ISOLATOR  IN  IMAGE  GUIDE 


ABSTRACT 

An  image  guide  consisting  of  a  layer  of  dielectric  and  a 
layer  of  ferrite  with  periodic  metal  strips  between  these  layers 
is  analyzed.  Non-reciprocity  of  the  phase  constants  caused  by  the 
gyrotropic  nature  of  the  ferrite  makes  it  possible  to  create  non¬ 
reciprocal  leakage  phenomena  by  means  of  the  metal  strips.  These 
phenomena  can  be  used  to  develop  a  distributed  isolator.  Numerical 
results  and  data  from  the  experimental  investigation  performed  are 
presented . 
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INTRODUCTION 


In  the  last  decade,  the  microwave  engineers  have  been 
trying  to  extend  integrated  circuits  at  microwave  frequencies  to 
millimeter-wave  frequencies.  The  approach  to  millimeter-wave  I.C.'s 
tends  more  to  optical  waveguides  than  to  coaxial  lines  as  is  the 
case  for  microwave  I.C.'s.  The  most  promising  means  to  make 
millimeter-wave  I.C.'s  is  the  image  guide  system.  Various  compo¬ 
nents  have  been  built.  However,  to  date  little  has  been  done  to 
create  a  non-reciprocal  device  compatible  with  image  guide  systems. 
The  structure  proposed  here  is  one  that  is  compatible  with  image 
guides  and  is  reasonable  in  size  at  millimeter-wave  frequencies. 

The  isolators  available  at  millimeter-wave  frequencies 
are  the  Faraday  rotation  type,  the  resonance  type  using  hexagonal 
ferrites,  and  those  made  from  a  junction  circulator.  None  of  these 
are  very  comp.ttible  with  image  guide  systems  although  the  resonance 
type  could  possibly  be.  Also,  hexagonal  ferrites  are  not  commer¬ 
cially  available. 

Several  structures  have  been  proposed  to  make  a  non- 
reciprocal  device  compatible  with  image  guides.  These  are  mainly 
field  displacement  devices  and  structures ^ ^ ^ that  have  physi¬ 
cal  resonances  at  different  frequencies  in  the  opposite  direc- 

f  1  g  1 

tions^  .  The  present  structure  is  an  extension  of  a  previously 
suggested  one  where  a  periodic  grating  was  created  in  the  dielectric 


1 
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guide  attached  to  a  ferrite  guide  .  The  structure  analyzed  in  this 
article  consists  mainly  of  a  ferrite  guide  next  to  a  dielectric 
guide  with  periodic  metal  strips  between  the  guides  and  oriented  in 
the  E-plane. 

In  this  article.  Chapters  1  and  2  review  ferrite  theory 
at  microwave  frequencies  and  with  a  simple  explanation  of  how  the 
non-reciprocal  ferrite  devices  presently  in  use  work.  Chapter  3 
deals  with  the  derivation  of  the  theory  for  analyzing  the  propaga¬ 
tion  constant  in  the  proposed  structure.  Chapter  4  presents  the 
numerical  procedure  and  the  computer  program.  Chapter  5  deals  with 
the  experimental  investigations  and  the  comparison  between  the 
experimental  and  the  theoretical  results. 


CHAPTER  1 


GENERAL  FERRITE  THEORY,  PARAMETERS  AND  MATERIALS 

The  purpose  of  this  chapter  is  to  give  some  insight  of 
what  ferrite  is,  what  makes  it  different  from  other  materials  and 
what  parameters  distinguish  different  types  of  ferrite  materials. 
This  chapter  also  deals  In  some  detail  with  the  derivation  of  the 
so  called  permeability  tensor  of  ferrites.  This  tensor  permeability 
is  the  most  direct  tool  which  could  give  an  engineer  insight  into 
how  ferrite  material  is  going  to  behave.  The  tensor  permeability 
is  what  makes  ferrites  different  from  other  dielectric  materials  and 
it  would  be  used  back  and  forth  in  analysis  of  systems  containing 
ferrite  materials,  so  it  would  be  a  good  idea  to  look  at  the  tensor 
and  see  how  it  originates  and  how  to  control  it  to  obtain  the 
desired  characteristics. 
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1.1  Types  of  Magnetism 

a.  Diamagnetism  and  Paramagentism 

For  any  material  the  magnetic  moment  of  an  atom  is  equal  to 
the  sum  of  the  moments  of  Its  electrons.  In  a  large  number  of  cases, 
this  moment  is  zero  due  to  the  symmetry  of  the  atom  which  contains  as 
many  electrons  rotating  in  one  sense  as  in  the  other  sense.  Just  as 
in  electromagnetism,  where  a  flux  gives  rise  to  an  induced  current 
tending  to  oppose  any  change,  a  substance  which  is  acted  on  by  a 
magnetic  field  tends  to  respond  by  the  creation  of  internal  charges. 
If  the  moments  of  the  atoms  are  zero,  then  the  induction  is  smaller 
than  the  applied  field  and  the  magnetic  susceptibility  must  be  nega¬ 
tive.  Such  a  substance  is  said  to  be  diamagnetic.  The  magnetic 
susceptibility  is  very  small  of  the  order  of  10  ^  (e.g.,  water  has 
X  =  -0.91  X  10"^) . 

If  the  magnetic  moment  of  the  atom  is  not  zero,  the  body 
is  said  to  be  paramagnetic.  Substances  where  atoms  have  an  odd 
number  of  electrons  arc  paramagnetic.  This  is  the  case  of  atoms 
which  have  an  incomplete  internal  shell  as  well  as  for  free  ions. 
Metals  are  often  paramagnetic  because  of  the  presence  of  conduction 
electrons,  as  are  also  some  bodies  which  have  an  even  number  of 
electrons. 

The  susceptibility  is  now  positive  but  it  is  still  very 
-5  -4 

small  (of  the  order  of  10  to  10  )  . 
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b.  Ferromagnetism 


A  substance  is  said  to  be  ferromagnetic  when  it  possesses 

a  spontaneous  magnetic  moment  in  the  absence  of  any  applied  magnetic 

field.  The  spontaneous  magnetic  moment  per  unit  volume  is  called 

the  saturation  intensity  M  . 

s 

An  explanatifin  of  this  phenomenon  was  first  given  by 

[29] 

VJeiss  who  postulated  the  existence  of  a  force  tending  to  align 
all  the  magnetic  moments  of  the  atoms  in  the  material  in  one  direc¬ 
tion.  Thus,  a  ferromagnetic  substance  is  a  paramagnetic  substance 
in  which  there  exists  a  field,  called  the  molecular  or  exchange 
field  by  Weiss.  This  field  aligns  the  magnetic  moments  of  all  the 

paramagnetic  atoms  into  one  direction. 

[22] 

Heisenberg  later  showed  that  the  molecular  field  was  a 
coupling  of  electrical  origin  which  existed  only  if  the  distances 
between  the  atoms  in  the  material  lay  between  certain  well-defined 
limits.  The  energy  of  this  coupling  is  called  the  exchange  energy. 
It  can  be  shown  that  for  two  neighboring  atoms,  whose  spin  moments 
equal  S  and  at  an  angle  <p  to  each  other,  the  exchange  energy  is 
given  by 


W  =  J  S 


2,2 


(1.1) 


where  J  is  a  coefficient  called  the  exchange  integral. 

However,  the  exchange  energy  which  tends  to  align  the 
moment  is  more  or  less  counterbalanced  by  the  thermal  energy  which, 
by  random  agitation  of  the  atoms,  tends  to  destroy  the  alignment. 
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This  energy  is  proportional  to  the  product  of  Boltzmann's  constant 

and  the  absolute  temperature.  At  low  temperatures  the  exchange 

energy  is  predominant  and  a  large  number  of  the  moments  are  aligned. 

The  saturation  intensity  M  is,  thus,  relatively  important.  As  the 

s 

temperature  rises,  the  increase  in  the  thermal  energy  gradually  des¬ 
troys  the  alignment  and  the  saturation  intensity  decreases.  When 
the  two  energies  are  equal,  the  saturation  intensity  becomes  statis¬ 
tically  zero.  The  temperature  at  which  this  occurs  is  called  the 
Curie  point  or  Curie  temperature. 

Above  the  Curie  point,  the  substance  is  no  longer  ferro¬ 
magnetic  and  becomes  paramagnetic.  Below  it,  the  moment  varies 
as  shown  in  Fig.  (1.1). 


Fig.  (1.1)  The  variation  of  saturation  magnetization  with  temperature 


The  value  of  the  molecular  field  can  be  deduced  from  a 
knowledge  of  the  Curie  temperature  T^.  In  the  case  of  iron  it  turns 


out  to  be  of  the  order  of  5  million  oersteds. 


The  exchange  energy  is  important  only  if  the  distances 


between  atoms  lie  between  very  narrow  limits.  There  exist  relatively 
few  ferromagnetic  materials.  In  fact,  only  iron,  cobalt  and  nickel 
exhibit  this  property  at  normal  temperatures. 

It  should  be  noted  that  a  ferromagnetic  material  Is  not 
necessarily  magnetized  to  saturation.  The  reason  is  that  in  the 
absence  of  an  external  magnetic  field,  the  substance  is  divided 
into  small  domains,  called  Weiss^^^^  domains  in  which  the  moments 
are  aligned  spontaneously  and  which  are,  therefore,  saturated.  But 
the  direction  of  magnetization  varies  from  one  domain  to  another 
so  that  the  overall  magnetization  can  be  quite  small  or  even  zero. 
Fig.  (1.2)  shows  schematically  the  distribution  of  Weiss  domains  in 
a  ferromagnetic  material  and  the  direction  of  magnetization  for 
each. 


Fig.  (1.2)  Weiss  domain  in  Ferromagnetic  materials 

This  microscopic  domain  structure  has  been  verified  photographi¬ 
cally.  It  is  assumed  that  each  domain  is  separated  from  its  neigh- 

[13 1 

bors  by  a  thin  region,  called  a  Bloch  wall,  in  which  the  mag¬ 
netization  changes  direction. 
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In  the  presence  of  an  external  field,  two  phenomena  can 
happen;  an  Increase  in  the  size  of  the  domains  in  which  the  magneti¬ 
zation  lies  in  the  direction  of  the  applied  field,  or  a  rotation  of 
the  magnetization  directions  in  the  various  domains  tending  to  align 
them  with  the  field.  The  first  of  these  two  effects  seems  to  exist 
for  relatively  weak  fields  and  the  second  for  fields  approaching 
saturation.  When  the  latter  condition  is  reached,  as  in  Fig.  (1.3), 
all  the  domains  are  aligned  with  the  applied  field. 


Fig.  (1.3) 

(I"*  Reversible  zone 

(2)  Increase  in  the  size  of  certain  domains 

(3)  Rotation  of  the  magnetization  of  the  domains 

If  the  applied  field  is  reduced  or  removed  after  a  speci¬ 
men  of  ferromagnetic  metal  has  been  magnetized  to  saturation,  the 
domains  do  not  revert  to  their  original  positions.  A  certain 
amount  of  magnetization  or  remanent  induction  B^,  is  left.  If  it  is 
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large,  we  have  a  permanent  magnet;  If  it  is  small,  the  material  is 
said  to  be  soft.  The  corrective  field  is  the  field  which  must  be 
applied  to  reduce  the  remanent  induction  to  zero. 

c.  Antiferromagnetism 

In  a  ferromagnetic  body,  magnetized  to  saturation,  all 
the  moments  of  the  atoms  are  oriented  parallel  to  the  applied  field 
when  the  temperature  is  much  below  the  Curie  point.  It  can  happen, 
however,  that  the  moments  form  an  antiparallel  arrangement,  as  shown 
in  Fig.  (1.4).  In  this  case  the  magnetization  remains  weak,  the 
susceptibility  being  then  of  the  same  order  of  magnitude  as  that  of 
a  paramagnetic  substance;  the  material  is  than  said  to  be  antiferro- 
raagnetic . 


J>0:  Ferromagnetic  J<0:  Antiferromagnetic 


Fig.  (1.4) 

It  can  be  shown  theoretically  that  ferromagnetism  corresponds  to  a 
positive  exchange  Integral  and  antiferromagnetism  to  a  negative 


exchange  Integral . 
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d.  Ferrlraagnetlsm 

Knowing  the  molecular  composition  of  a  material,  one  can 
determine  the  number  of  magnetons  available  to  the  part  in  the  mag- 

[131 

netization;  (one  magneton  usually  defined  as  the  Bohr^  ■*  magneton 

is  the  smallest  magnetic  moment  possible,  which  occurs  due  to  one 

electron  spin) .  It  is  then  easy  to  deduce  the  saturation  intensity 

of  the  material  when  all  the  magnetons  are  aligned  parallel  to  each 

other.  It  is  now  found  that  for  certain  compounds  such  as  iron 

oxide,  FeO,  Fe20,  for  example,  the  saturation  intensity  obtained 

experimentally  is  considered  less  than  the  predicted  value.  If  the 

calculation  is  done  in  reverse,  it  is  found  that  only  4.2  magnetons 

per  molecule  of  liave  contributed  to  the  magnetization,  as 

against  the  14  magnetons  available.  To  explain  this,  Neel^^^^  pre- 

3+ 

sumed  that  the  trivalent  Fe  ions  of  Fe20^  were  antiferromagnetic 

2+ 

(moment  aligned  antiparallel),  whereas,  the  divalent  Fe  ions  of 
FeO  were  ferromagnetic  (moment  parallel) ;  analogous  ideas  had  been 
proposed  by  Gullland  on  the  subject  of  the  components  of  manganese 
and  bismuth.  This  type  of  magnetism,  which  applies  to  the  ferrites 
which  we  are  about  to  discuss,  has  been  given  the  name  ferrlmagnet- 
ism. 

In  a  ferrlmagnetlc  molecule,  then,  there  are  atoms  between 
which  the  exchange  integrals  are  positive  and  atoms  between  which 
they  are  negative.  The  overall  moment  of  the  molecule  is  the  dif¬ 
ference  between  the  sums  of  the  moments  of  opposite  senses. 


'.r<- 


1.2  Material  Composition,  Saturation  and  Curie  Temperature 

Relationships 

As  the  temperature  of  ferrite  is  raised,  the  magnetic 
structure  is  weakened  and  eventually  falls  into  complete  disarray 
(the  Curie  point) .  When  non-magnetlc  ions  are  substituted,  the 
structure  is  already  weakened  and  the  Curie  point  reduced  Fig. 

f  2Q1 

(1.5)  .  Rare  earth  substitutions  for  yttrium  in  YIG  do  not 

weaken  the  interaction,  and  the  saturation  magnetization  can  be 
varied  without  affecting  the  Curie  temperature. 


Fig.  (1.5)  Saturation  Magnetization  and  Curie 

temperature  for  various  ferrite  systems 
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For  ferrite  devices,  the  Curie  point  determines  both  the 
maximum  operating  temperature  and  a  critical  temperature  range  where 
small  temperature  variations  will  have  large  effects  on  device 
characteristics. 

1.3  Precession  and  Gyromagnetic  Resonance 

Consider  a  particle  rotating  about  itself  (in  this  case 
a  spinning  electron)  which  we  shall  assume  to  have  the  shape  of  a 
small  sphere  Fig.  (1.6).  The  axis  of  rotation  is  along  a  direction 
Otj^,  which  makes  an  angle  0  with  the  Oz  axis  of  a  set  of  reference 
axes  with  origin  0  at  the  center  of  the  electron.  Since  the  elec¬ 
tron  has  mass,  there  exists  along  the  Ot  a  kinetic  moment  which  is 
designated  by  the  vector  p. 


Fig.  (1.6) 


at,**. 


It  can  be  shown  that  the  kinetic  moment  is  the  product  of 

the  action  quantum  h  and  the  spin  quantum  number  s,  which  has  only 

—  26 

two  possible  values  ±  1/2.  Since  h  =  105.4  x  10  J/s,  the  only 

two  possible  values  of  the  kinetic  moments  are  p^  =  hs  = 

—  36 

±  52.7  X  10  j/s.  We  know,  however,  that  the  electron  also  pos¬ 
sesses  a  magnetic  moment  M  which,  in  this  particular  case,  is  the 

-21 

Bohn  magneton  M^=9.27xlO  Mx-cm. 

Any  particle  having  mass  and  spin  will  give  simultaneously 
a  kinetic,  moment  and  a  magnetic  moment.  By  definition,  the  ratio 
Y  =  M/P  is  called  the  gyromagnetlc  ratio  of  the  particle.  In  the 
particular  case  of  the  electron: 


|y(  =  ^  =  17.6  X  10^ 
e 


If  no  external  force  is  applied  to  the  particle,  it 
rotates  about  the  axis  Ot,  where  direction  remains  fixed.  Let  us 
suppose  that  a  d.c.  magnetic  field  H  is  applied  along  the  Oz  axis. 
This  magnetic  field  exerts  on  the  moment  M  a  couple: 


c  =  M  X  H, 


(1.3) 


and  the  time  derivative  of  the  kinetic  moment  is: 


=  M  X  H 
dt 


(1.4) 


But  from  Y  =  M/p  we  also  get 


iE  =  i 

dt  Yd 
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t 


I 

I 

i 


Combining  the  above  equations  we  get 


1  ^ 
Y  dt 


=  M  X  H 


(l.G) 


which  is  t  iG  differential  equation  of  motion  for  the  moment  M. 

With  M,  the  tip  of  the  vector  M,  as  origin,  let  us  take 
a  set  of  coordinate  axis  Mt^^,  Mt2,  Mt^:  Mt^  is  in  the  place  of  M 
and  02,  and  Mt^  is  perpendicular  to  this  plane.  The  projections 
of  dPl/dt  on  these  three  axes  are  given  by  kinematics.  They  are 


on  Otj^ 


dt 


on  Ot^ 


MdO 

dt 


on  Ot^ 


M  sln6 


d(i> 
dt  ’ 


(1.7) 


being  the  angle  between  the  projection  of  M  on  xOy  and  the  Ox 
axis. 


Furthe 

three  axes  are: 

the  equations  of 

MH  ,  s  in  0 .  Th  1  s 
dc 


r,  the  projections  of  the  vector  product  on  the 

on  Ot,  &  Ot_  =  0;  and  on  0t_  =  MH,  sin9.  Finally, 
12  J  dc 

1  dM  1  „  d(J) 

mot  i(>n  of  M  are  written  —  —  =  0,  —  M  sin0  —  = 

Y  dt  Y  ut 

last  relationship  gives: 


d4> 

dt 


yH 


dc  ‘ 


(1.3) 


Since,  in  addition,  =  0  and  M  =  0,  then  the  locus  of  the  tip 

at  at 

of  the  M  vector  is  a  circle  lying  in  a  plane  perpendicular  to  Oz. 
The  circle  is  described  with  an  angular  velocity: 


15 


_ 

dt  ' 


(1.9) 


Tlie  M  vector,  therefore  describes  a  cone  around  the  Oz 
axis  with  semi  apex  angle  0.  We  recognize  here  the  phenomenon  of 
precession  in  a  gyroscope.  The  frequency  of  rotation  is  given  by: 

Y 

w  =  27Tf  SO  that  f  =  —  H.  .  In  the  case  of  an  electron, 

2ti  dc  ’ 

f  =  2.8H,  (1.10) 

dc 

where  f  is  in  GHz  and  H  is  kilo  oersteds. 

The  precession  thus  produced  is  accompanied  by  a  dissipa¬ 
tion  of  energy  which  appears  almost  entirely  as  heat,  the  remainder 
being  transformed  into  an  electromagnetic  radiation  of  frequency  f. 
This  dissipation  of  energy  brings  about  a  decay  of  the  motion,  the 
angle  6  becomes  progressively  smaller  and  precession  ceases  when  M 
is  aligned  with  H^^ .  The  time  between  the  application  of  and 
the  end  of  the  precessional  motion  is  called  the  relaxation  time. 

Let  us  now  consider  an  electron  under  the  influence  of  a 
d.c.  magnetic  field  for  a  time  long  enough  for  the  moment  M  to 
be  aligned  with  H  and  for  all  precession  motion  to  have  ceased. 

Fig.  (1.7). 


Fig.  (1.7) 


Let  us  apply  to  the  electron  a  second  magnetic  field  h, 
lying  In  the  xOy  plane.  The  resultant  field  will  be  a  vector 
+  h  making  an  angle  6  with  the  direction  Cz  of  tlie  moment  M, 
and  we  find  ourselves  in  the  situation  of  the  preceding  paragraph 
and  the  spin  moment  will  again  acquire  a  precessional  motion  around 
the  direction  of  (M  +  h) . 

Let  us  now  suppose  that  the  field  h  is  a  high  frequency, 

circularly  polarized  field.  It  is  known  that  such  a  field  is 

characterized  by  a  vector  h  which  rotates  in  the  plane  of  magnetic 

polarization  (here,  the  plane  xOy)  at  the  given  frequency.  In 

these  conditions,  the  vector  (H,  +  ii)  describes  a  cone  around  Oz , 

dc  ’ 

carrying  with  It  the  moment  M  which  acquires  a  precessional  motion 
around  it.  Two  cases  can  then  arise: 
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a.  The  sense  of  rotation  of  h  Is  such  that 

the  cone  is  described  by  (H ,  +  h)  in 

dc 

the  same  sense  as  that  of  the  precession; 
in  this  case,  if  the  frequency  of  h  Is 
the  same  as  that  of  the  precession,  there 
will  be  synchronization  of  the  two  phenomena 
and  we  have  resonance.  The  amplitude  of 
the  precession  (aperture  of  the  cone) 
Increases  until  the  energy  supplied  by 
the  high  frequency  field  exactly  balances 
the  losses,  and  the  motion  is  maintained 
in  equilibrium  at  the  expense  of  the  high 
frequency  field  which,  therefore,  suffers 
a  significant  absorption.  This  phenomenon 
is  known  as  gyromagnetic  resonance. 


If  h  is  small  compared  with  H,  it  appears 
when  the  value  of  the  latter  is  such  that 
the  frequency  given  by  f  =  2.8H^^  is  equal 
to  the  frequency  of  the  field  h.  The  d.c. 
field  at  resonance.  Ho,  is  therefore 


iio 


f 

2.8 


b.  The  sense  of  rotation  of  h  is  such  that 
the  cone  is  described  by  +  R)  in 

the  opposite  sense  to  that  of  the  precession. 
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In  this  case,  synchronization  is 
impossible  and  there  cannot  be  any 
gyromagnetic  resonance. 

One  important  fact  will  be  observed:  resonance  occurs 
for  only  one  sense  of  polarization.  Now,  for  a  circularly  polar¬ 
ized  high  frequency  wave,  the  sense  of  rotation  depends  on  the 
direction  of  propagation.  Resonance,  therefore,  occurs  in  one 
direction  of  propagation  only;  in  the  opposite  direction,  no  signi¬ 
ficant  absorption  occurs.  This  is  a  formal  contradiction  of  the 
reciprocity  theorem.  Gyromagnetic  resonance  is,  therefore,  a  non- 
reciprocal  phenomenon. 

The  absorption  curve  of  ferrite  material  vs.  will  be 

as  shown  in  Fig.  (1.8).  One  may  notice  that  changing  the  direc¬ 
tion  of  H  ,  or  the  microwave  field  has  just  the  same  effect, 
dc 


Fig.  (1.8) 
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l.A  The  Permeability  Tensor  and  Line  Width 


In  the  simplest  microwave  case,  the  total  effective 

dii  - 

magnetic  field  in  the  equation  of  motion  (  —  =  Y  M  x  H  )  consists 
of  the  d.c.  magnetic  field  and  the  r.f.  magnetic  field  h. 


H  =  +  h 

dc 


(1.11) 


The  total  magnetization  consists  of  the  d.c.  magnetization  ^ 
and  the  r.f,  magnetization  m 


M  =  M,  +  m  . 

d.c. 


(1.12) 


In  component  form,  the  above  equations  are: 


dc 


r  *1 

. 

- 

0 

h 

0 

in 

X 

X 

0 

,  h  = 

h 

y 

’  ”dc  = 

0 

,  in  = 

m 

y 

H 

h 

m 

dc 

z 

dc 

z 

- 

- 

-  ^ 

The  equation  of  motion  now  gives 


dm 

=  m  Y  (H,  +  h  )  -  h  Y  (M,  +  m,) 

dt  y  dc  z  z  dc  z 


dm 

=-m  y(H,  +h)+h  y(M,  +m) 
dt  X  dc  z  X  dc  z 


(1.13) 


dm 

-r—  =  m  Yh  -  m  Yh 
dt  X  y  y  X 


■I 
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In  the  small  signal  approximation  higher  order  terras  of 
m  and  h  are  set  equal  to  zero.  The  small  signal  approximation  is, 
therefore; 


dm 

-7-^  =  m  yH  -h  yM, 
dt  y  dc  y  dc 


dm 


— f  =  -m  Y  H  +  h  Y  M , 
dt  X  dc  X  dc 


(l.U) 


dm 

z 


=  0 


have: 


Rewriting  and  differentiating  the  last  equations,  we 


m+w  m=wuh-wh 
X  ox  mox  my 


m+w  ra=ioh+u)wh 
y  oy  mx  moy 


(1.15) 


m  =  0  , 

z 


where  oj  =  yM  ,  ,  w  =  yH  ,  ,  and  •  indicates  time  derivative, 

m  dc  o  dc 

If  the  time  dependence  of  the  r.f.  quantities  is  of  the 
form  exp(jojt),  a  susceptibility  tensor  [y]  can  be  defined  which 
relates  the  r.f.  magnetization  to  the  r.f.  magnetic  field. 


m  =  Ixl  h 


Where  [xl 


'S.X  \y  0 

^X  Xyy  0 

0  0  0 


(1.16) 


-V"  ^ 


r 
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and  X 


XX 


U  (U 

m  o 

2  2 
0)  -  to 

o 


1(0  (0 
m 

2  2 
to  -  to 
o 


(1.17) 


(1.18) 


The  components  of  the  susceptibility  tensor  have  a 
singularity  at  to  =  to^.  This  is  defined  as  the  resonance  condition, 
A  tensor  permeability  can  also  be  defined  by  relating 
the  r.f.  flux  density  b  to  the  r.f.  magnetic  field  ii, 


b  =  y  h  +  m 
o 


(1.19) 


or 


b  =  u  [y^]  ii 

o  r 


where  [y^]  =  [1]  +  [ xJ • 

In  component  form,  the  tensor  permeability  is  given  by 


[yrl 


■  y  jk  O' 
jk  y  0 
.0  0  1  . 


where  y  =  1  + 


(1.20) 


jk  =  -x^y 

[21] 

This  is  the  well  known  Polder  tensor  permeability. 
The  permeability  tensor  is  usually  employed  in  conjunction  with 
Maxwell's  equations  to  solve  microwave  ferrite  boundary  problems. 


Damping  may  be  introduced  into  the  equations  of  y  and  k  by  adding 


an  imaginary  frequency  term  jam  to  in  these  equations.  We  shall 
next  discuss  where  the  jam  term  comes  from  and  what  is  the  value  of 

a. 

To  stabilize  the  motion  of  the  magnetization  vector  at 
resonance,  a  damping  term  must  be  Introduced  into  the  equation  of 
motion.  One  form  of  phenomenological  damping  term  is  due  to  Landau 
Lifshitz^^^^  given  by: 


^  =  y(M  X  H)  -  [M  X  (M  X  H)]  (1.21) 

where  a  is  a  dimensionless  parameter  chat  determines  the  damping. 
Solving  for  y(M  x  H)  in  the  above  equation  gives 

•>(M  X  H)  =  ~  +  iM  X  (M  X  ii)  ]  (1.22) 

and  introducing  this  into  the  last  part  of  the  previous  equation 
gives ; 


dt 


y(M  X  H) 


1m[  ^  ^  ■ 


ya 


M  X  [M  X  (ii  X  H)  ]  .  (1.23) 


For  most  ferrites  the  last  term  is  negligible  since  a  is  small. 
Therefore,  the  equation  could  be  written  as, 

-^(Mxf).  (1.2M 

This  damping  term  is  shown  in  Fig.  (1.9).  The  components  of  the 
tensor  susceptibility  obtained  from  the  small  signal  solution  to 
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the  above  equation  containing  the  damping  terras  are: 


0.'  (a  +  ioia) 
mo 


(1.25) 


1  0)  U) 

m 


\x  \ 2  2 

(to  +  jtoa)  -  (1) 


(1.26) 


Hence,  the  only  difference  between  these  equations  and  the  equations 

when  we  did  not  consider  damping  is  that  an  imaginary  frequency  term 

joja  has  been  added  to  the  resonant  frequency  to^.  Hence,  the  damping 

term  can  always  be  introduced  by  replacing  w  by  ui  +juu  in  the  loss  free 

components.  The  real  and  imaginary  parts  of  the  susceptibility 

tensor  are  plotted  in  Fig.  (1.10)  as  a  function  of  lo  /lo,  u  /lo  and  a. 

mo 
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Thf  damping  factor  a  can  be  related  to  the  linewldth  of 
ferrites  in  the  following  way.  The  linewidth  is  usually  defined 
as  the  difference  between  the  magnetic  field  values  at  a  constant 
frequency  where  the  Imaginary  parts  of  the  diagonal  component 

of  the  susceptibility  tensor,  attains  the  value  of  half  its  value 
at  resonance  as  shown  in  Fig.  (1.11). 


Fig-  (1.11)  Linewidth  of  Ferrites 


Kquating  the  imaginary  part  of  to  half  the  maximum  value  of 

■j  -> 

(u  ) 


2  2 

W  Wft  (w  ,  -  0)  ) 

m  Ui 


2  2  2  2  2  2 
-  u)  )^  -F  a 


=  1/2 


2ioa 


(1.27) 


Rearranging  gives  a  quadratic  in  w 


01 


01 


2  2  4  2 

-  2(1)  (Dqj^  +  (d  (1  -  4a  )  =  0, 


4 


(1.28) 
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whose  solution  Is 


Hence , 


1/2 

2  “  '*'^1  -  2a)  =  w  ±  ua. 


“02  “  “oi 


(1.29) 


(1.30) 


The  linewidth  is  usually  defined  by 


AH  - 


“02  “Ol 


(1.31) 


yAH 

Uu)  =  . 


(1.32) 


For  engineering  purposes  the  imaginary  frequency  jaio  is  usually 
written  in  terms  of  AH  from  the  last  equation. 


1.5  The  Permeability  Tensor  Below  Saturation 


Below  saturation  there  exists  no  additional  formulation 
of  the  tensor,  since  not  all  the  electron  spins  are  aligned  to  the 
applied  magnetic  field,  and  we  do  not  have  any  measure  of  alignment 
to  formulate  a  theory  which  would  give  us  a  good  result. 

The  only  remaining  solution  is  to  measure  the  tensor  and 
then  fit  empirical  equations  to  the  experimental  results.  One  such 
work  has  been  done  by  J.  J.  Green  and  F.  Sandy^^^^.  In  their  1974 
paper  they  discuss  the  experimental  results  of  a  number  of  different 
types  of  ferrites  and  then  fit  the  appropriate  general  empirical 
formulas  which  are  given  by 


(1.33) 


where 


U-  =  2/3  [1  -  (to  +  1/3 

u  m 


p  =  Uq  +  (1  -  Pp)  (0)^/0)^' 


(1  -  w  /o)  )■ 

p^  =  Py  0  m 


Measurement  of  the  Permeability  Tensor 


The  permeability  tensor  obtained  by  measurement  is  perhaps 
the  best  way  to  characterize  a  ferrite  material.  Actually  when 
ferrites  are  produced  the  only  way  to  obtain  important  parameters 
such  as  saturation  magnetization  and  linewidth  is  by  making  mea¬ 


surements. 


There  exist  various  types  of  measurement  techniques,  the 
main  two  being  the  cavity  perturbation  and  the  ferrite  filled  wave¬ 
guide  method.  The  former  method  is  usually  preferred  since  it 
gives  better  results  and  the  sample  size  required  is  very  small. 

The  cavity  perturbation  technique  could  be  classified  into  two 
methods,  namely  the  non-degenerate  cavity  method  and  the  degenerate 
cavity  method  due  to  J.  0.  Artman^^^.  Using  non-degenerate  cavities 
provides  information  only  about  the  diagonal  element  of  the  tensor. 


However,  at  times  this  is  sufficient,  especially  close  to  resonance 
and  for  small  damping  where  the  diagonal  and  off-diagonal  elements 
are  closely  related  by  equation 


28 


k  =  1  -  p 


(1.34) 


Degenerate  cavity  methods  give  Information  about  the 
diagonal  and  off-diagonal  element  of  the  tensor  but  the  overall 
inaccuracies  involved  in  having  a  good  degenerate  cavity  make  it 
difficult  to  obtain  reasonable  results  with  this  method. 

The  perburbation  theory  for  various  types  of  cavities 
and  S'’"  pie  shapes  have  been  derived;  a  good  account  of  this  is  given 
by  Lax  and  Button^ A  general  idea  of  the  derivation  of  per¬ 
turbation  relations  is  given  next. 

In  the  cavity  methods  the  way  to  determine  the  suscepti¬ 
bility  tensor  Is  to  measure  the  shift  in  resonant  frequency  and 
the  change  in  Q  of  the  cavity  loaded  by  a  small  sample  of  ferrite 
material . 


In  the  case  of  the  unperturbated  cavity  resonating  at  a 
frequency  (o^  Maxwell's  equations  are 


7  X  Eq  =  -ju,p  Pp  hp 
V  X  ho  =  ja,o  to  % 


(1.35) 

(1.36) 


If  the  cavity  is  now  perturbed  by  a  sample  having  a 


volume  Av,  the  cavity  will  resonate  at  a  new  resonant  frequency  o). 
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Maxwell's  equations  applicable  to  the  perturbed  cavity  are 

A  X  E  =  -jwiiQh  (1  +  [xD  (1.37) 

A  X  h  =  joseQE  (1  +  [Xj^]).  (1.38) 

[  x]  and  [Xj^]  are  in  general  tensors,  although  in  a  ferrite  [Xj^]  is 
a  scalar  Xj^. 

In  the  above  equations  the  electric  and  magnetic  suscepti¬ 
bility  terms  are  zero  everywhere  outside  the  region  Av  occupied  by 
the  ferrite  material.  After  some  algebraic  manipulations,  and 
taking  into  account  the  boundary  conditions  at  the  surface  of  the 
cavity,  we  obtain 


CO  -  ai_ 


“0 


~-l/v  ‘^0  f 


(1.39) 


This  equation  is  exact  and  could  be  evaluated  if  the  field  compo¬ 
nents  of  the  perturbed  cavity  were  known.  However,  according  to 
first  order  perturbation  theory,  we  can  set  E  ==  E^,  and  h  =  h^ 
tliroughout  the  cavity. 

If  the  sample  is  located  at  a  point  in  the  cavity  where 
Eq  =  0,  then  we  have 


-j/v  (iXtl-^j-hpldv 

2  /v  "o  • 


J  dv 


(1.40) 
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in  tlie  above  equation  we  have  used  the  fact  that  the 
stored  electric  and  magnetic  energies  averaged  out  over  time  are 
nearly  equal.  We  have  also  made  use  of  the  fact  that 


lx]  ii  =  [x^]  hp 

This  equation,  therefore,  gives  the  external  susceptibility  tensor 
of  the  material.  Since  the  susceptibility  is  complex  in  general, 
the  frequency  oj  is  complex  also.  For  small  change  this  can  be 
written  as 


where  is  the  unloaded  Q  of  the  cavity.  Hence, 

%  .  -X  ‘'"d  •  V  •  *'0  ““ 

X  ■  'o> 

,,  1 ,  Xv  <i«d"  -"X" 

y.  %  ‘  V 


(1.41) 


(1.42) 


(1.43) 


The  real  and  imaginary  parts  of  the  susceptibility  tensor  are, 

therefore,  obtained  by  measuring  the  shift  in  resonant  frequency 

and  the  change  in  unloaded  Q  of  the  cavity.  A  detailed  description 

[24] 

of  how  to  perform  such  measurements  is  given  by  Sucher  and  Fox 


— r- 
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1.7  Commercially  Available  Ferrites 


Properties  of  typical  commercially  available  fen ite 
materials  are  shown  in  Figs.  (1.12)  and  (1.13)^^^’^^’^.  The  satura¬ 
tion  magnetization  varies  from  250G  for  substituted  YIC  to  5000G 
for  NiZn  ferrites;  Mg  Mn  A1  ferrites  have  a  saturation  which  ranges 
from  600  to  3000G.  Actually  these  are  the  general  purpose  ferrites 
mainly  used  at  X-band. 

Materials  with  narrow  linewidth  will  be  mainly  substituted 
YTG  or  Mg  Mn  A1  ferrites.  Addition  of  Ni  or  Co  usually  tends  to 
increase  the  linewidth.  The  dielectric  loss  tangent  tends  to  be 
lower  for  garnets  and  MG  ferrites  than  for  Mn  and  Ni  ferrites  so 
if  a  high  saturation  is  not  required,  it  appears  to  be  advantateous 
to  choose  garnets  or  Mg  Al  ferrites  for  low  magnetic  and  dielectric 
loss  over  a  wide  range  saturation  magnetization. 
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Fig.  (1.12) 


CHAPTER  2 


CONVENTIONAL  NON-RECIPROCAL  EFFECTS 

Before  describing  the  original  work  done  in  this  research 
project,  it  is  worth  mentioning  what  has  already  been  developed  in 
the  field  of  non-reciprocal  devices.  A  brief  description  of  most 
of  the  non-reciprocal  effects  discovered  up  to  now  is  given  in 
this  cahpter.  Finally,  a  brief  discussion  on  the  merit  and  improve¬ 
ment  methods  of  each  type  of  device  will  be  given. 
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2.1 


Resonance  Effect 


Consider  a  piece  of  ferrite  in  which  a  circularly  polar¬ 
ized  wave  is  propagating  in  a  given  direction.  Let  us  suppose  that 
the  plane  containing  tlie  magnetic  field  vector  of  the  wave  is  per¬ 
pendicular  to  a  d.c.  magnetic  field  applied  to  the  fe'rrite.  These 
are  the  conditions  under  which  resonance  is  possible  and,  if  the 
circular  polarization  is  rotating  in  the  positive  direction,  we 
obtain  the  maximum  value  of  ()i"+)  at  a  certain  value  of  the  d.c. 
applied  field.  For  a  wave  propagating  in  the  leverse  direction,  it 
is  (p"-)  which  is  Involved  and  we  know  that,  under  these  conditions, 
quite  a  large  attenuation  is  possible  in  the  ferrite. 

It  is  easy  to  visualize  an  application  of  this  phenomenon 
with  the  arrangement  shown  in  !'ig.  (2,1).  In  a  circular  guide,  a 
wave  is  propagating  with  circular  polarization  in  the  mode. 


Forwird  direction  Reverse  direction 

Fig.  (2.1)  Circular  guide  filled  with  ferrite 
and  magnetized  along  the  axis 


■niiwiy-yy. 
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In  a  rectangular  guide,  when  propagation  is  in  the  mode,  the 

components  of  the  magnetic  field  in  the  directions  Ox  and  Oy  are  as 
follows; 

h  =  sin^e-^^y  (2.1) 

X  L 


where 


h 

y 


n 

L 


cos 


(2.2) 


The  amplitudes  of  these  components  depend  only  on  x,  while 
the  exponential  term  is  representing  the  phase  of  the  propagating 
field  component.  Since  circular  polarization  is  defined  by  two 
components  in  quadrature  in  space  and  90°  apart  in  phase,  for  all 
points  where; 


li=jh  (2.3) 

X  y 

the  polarization  of  the  magnetic  field  is  circular. 

This  relationship  is  verified  for  a  plane  parallel  to 
the  narrow  side  of  the  waveguide  at  a  distance  x^  from  the  origin 
such  that: 


-  j  6  s  in 


itxo 


.  TT  Tnto 

_j  _  eos  — 


tan 


TOO  Ag 

'  1.  °  2L 


But  in  a  rectangular  guide: 

Xg  = 


Xo 


/U  -  (Xo/2L)2] 


(2.4) 

(2.5) 


(2.6) 


So  the  formula  shows  that  xo  depends  on  frequency.  Using  this 
formulation,  we  could  locate  the  position  where  we  have  circular 
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polarization  in  a  rectangular  waveguide.  A  piece  of  ferrite  will  be 
placed  at  this  position  in  the  waveguide  such  that  the  effect  same 
as  in  the  circular  waveguide  will  occur.  For  an  applied  static 
magnetic  field,  the  ferrite  material  will  go  into  resonance  and 
ahsorb  the  microwave  energy  for  one  direction  of  propagation  while 
nothing  will  happen  in  the  opposite  direction.  The  main  problem  of 
such  a  device  is  that  the  position  where  we  have  circular  polariza¬ 
tion  varies  with  frequency.  Various  techniques  are  used  to  counter¬ 
act  this  shift  with  frequency.  The  most  popular  methods  used  are 
placing  a  piece  of  dielectric  next  to  the  ferrite  or  by  tapering  the 
ferrite  cross  section. 

2.2  The  Faraday  Effect 

Consider  again  the  arrangement  in  Fig.  (2.1),  but  suppose 
now  that  a  TE^^^  wave  with  rectilinear  polarization  is  propagating 
along  the  guide  from  left  to  right. 

We  know  that  a  linearly  polarized  wave  can  be  resolved 
into  two  circularly  polarized  waves  of  the  same  amplitude  and  oppo¬ 
site  sense  of  rotation.  In  the  ferrite  acted  on  by  the  field 
one  of  the  two  waves  will,  therefore,  be  affected  by  the  permeabil¬ 
ities  p+  and  the  other  by  the  permeabilities  p-.  The  resolution 
of  the  linear  polarization  into  two  circular  polarizations  is 
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represented  in  Kig.  (2.3(a)).  It  can  be  seen  that  the  vector  R  is 
the  sum  of  the  two  vectors  and  rotating  in  opposite  directions. 


Fig.  (2.3)  The  vector  representation  of  the  Faraday 
effect,  considering  the  wave  propagating 
down  the  guide  as  two  circularly  polarized 
waves  in  opposite  directions 


and  being  synchronous,  R  remains  always  vertical. 

If,  in  the  course  of  propagation,  one  of  the  two  components  slowly 
changes  phase  with  respect  to  the  other,  R  will  rotate  as  shown  in 
Fig.  (2.3(b)).  The  component  P^ ,  instead  of  being  at  P^ ,  symmetri¬ 
cal  about  the  vertical  with  P^,  lags  by  an  angle  i}!  so  that  the  posi¬ 
tion  of  R  is  now  as  indicated  on  the  diagram  and  makes  an  angle  6 
with  the  vertical. 

Now  this  is  exactly  what  happens  in  the  guide.  Since  p'+ 
and  p'-  are  different,  the  propagation  constants, 

2ll  j  D~  2tT  /o  7'! 

B  =  — :j:  and  6  =  — ^  (2.7) 

are  different. 


-5Tr 
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The  relative  lag  of  one  component  behind  tlie  other  is; 


(t>  =  t3‘*'  -  e]  t 


(2.8) 


so  that,  the  rotation  of  the  direction  of  R  is: 

e  .  ■ 


(2.9) 


This  rotation  of  the  plane  of  polarization  of  the  wave, 
which  is  also  found  in  optics  is  the  Faraday  effect.  Faraday  rota¬ 
tion  depends  on  the  relative  values  of  p'"*”  and  p'  and,  to  make  use 
of  it,  it  is  necessary  to  work  in  a  region  where  they  differ  appre¬ 
ciably.  However,  ti"^  and  p"  must  be  small  so  that  the  system  will 
not  introduce  any  appreciable  attenuation  in  the  propagation. 

Tlie  essential  feature  of  the  Faraday  effect  is  its  non- 
reprocity.  Let  us  assume  that,  in  the  case  of  Fig.  (2.1),  for  prop¬ 
agation  from  left  to  right,  the  Faraday  rotation  takes  place  from 
left  to  right  as  seen  by  an  observer  looking  in  the  direction  of 
propagation  of  the  wave.  In  the  reverse  direction  of  propagation, 

the  position  would  be  the  same  if  H,  were  also  reversed;  but  H, 

dc  dc 

has  maintained  its  direction  so  that  the  p  are  permuted  and  Faraday 
rotation  takes  place  from  right  to  left  as  seen  by  an  observer 
looking  in  the  direction  of  the  Inverse  wave  and,  therefore,  from 
left  to  right  as  seen  by  the  first  observer.  The  two  rotations 
corresponding  to  the  forward  and  reverse  directions  are  equal  and 
add;  it  suffices,  therefore,  to  make  0  =  45®  to  produce  an  isolator. 


AO 


The  structure  used  is  that  of  Fig,  (2. A)  which  shows  also 
the  position  of  the  electric  field  vectors.  A  rectangular  guide 
excites  a  circular  guide  containing  the  ferrite.  The  axial  magne¬ 
tic  field  is  created  by  a  coil  around  the  circular  guide.  The  lat¬ 
ter  in  turn  excites  a  rectangular  guide  set  at  AS®  to  the  input 
guide.  It  can  be  seen  that  for  a  Faraday  rotation  of  A5°  there  is 
transmission  from  left  to  right.  In  the  r Ight-to-lef t  direction, 
the  rotation  brings  the  electric  field  into  quadrature  with  the 
field  of  the  TE^^  mode  in  the  input  guide.  The  field  in  quadrature 
cannot  be  propagated  in  the  guide.  It  is,  therefore,  reflected  to 
the  right  where  after  suffering  a  rotation  of  A5“  it  finds  itself  in 
quadrature  with  the  normal  field  of  the  output  guide  from  which  it 
is  again  reflected  to  the  left  and  again  rotated,  this  time  into 
the  right  direction  to  enable  it  to  re-enter  the  input  guide.  To 
prevent  this  an  absorbent  load  is  placed  near  the  input  guide;  this 
load  is  a  resistive  slab  perpendicular  to  the  input  field  (and  so, 
having  no  effect  on  it)  and,  consequently,  parallel  to  the  field 
coming  from  the  right,  thus  insuring  maximum  absorption. 
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Fig.  (2.4)  The  Faraday  rotation  isolator. 

a .  Cutaway  view 

b.  Method  of  operation 


It  will  be  seen  that  the  ferrite  employed  is  a  cylindrical 
rod  which  occupies  only  a  small  part  of  the  cross-section  of  the 
circular  guide.  In  fact,  a  small  cross-section  of  ferrite  is  suf¬ 
ficient  to  affect  a  large  part  of  the  transmitted  energy  because  of 
the  high  values  of  e  and  p;  it  Is  found  that  there  is  an  optimum 
ratio  of  the  diameters  of  the  ferrite  and  the  guide  for  a  maximum 


figure  of  merit. 
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2.3  Ttu'  Field  Displacement  Fffect 

Tlie  field  displacement  effect  Is  due  to  the  appearance  of 
liic.lur  modes  in  the  ferrite  as  a  result  of  the  liigii  value  of  tlie 
product  of  its  dielectric  constant  and  its  permeability.  Since  the 
pi  rmeab il  ity  tensor  depends  on  the  direction  of  propagation  of  the 
wave,  it  is  conceivable  that  one  might  get  non-reciprocal  effects. 

Tlie  mechanism  of  field  displacement  can  he  represented  by 
Fig.  (2.'))  for  a  rectangular  guide  loaded  with  a  slab  of  ferrite 
stuck  to  oni'  of  the  narrow  sides  of  the  gutdi'. 


Fig.  (2.5)  Kieliis  in  a  wave  guide 

partially  filled  with  ferrite. 

a.  The  field  distribution  of  the 
two  modes 

b.  Resultant  field  in  +ve  direction 

c.  Resultant  field  in  -ve  direction. 


In  this  guide,  partially  filled  with  dielectric  (ferrite),  there  Is 
a  mode  whose  configuration  is  F. :  this  mode  is  Independent  of 


the  direction  of  propagation  and  is  a  function  of  the  dimensions  of 
the  ferrite  and  its  dielectric  constant  only.  However,  the  high 
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value  of  the  permeability  tensor  causes  higher  modes  to  appear,  the 

effect  of  which  is  indicated  by  E^,  the  values  of  E^  depending  on 

the  direction  of  propagation.  If  we  denote  by  E^  the  distribution 

n 

obtained  for  one  direction  and  by  E^  that  obtained  for  the  opposite 

direction,  the  resultant  field  is  E  +  e"*"  in  the  one  case  and  E  +  E 

n  n 

in  the  other. 

The  resultant  configurations  are  arranged  as  shown  in 
Fig.  (2.5(b))  and  (2.5(c)).  In  the  former  case  nearly  all  the 
energy  is  concentrated  in  the  ferrite  with  a  field  maximum  at  the 
air-ferrite  Interface,  whereas,  in  (c)  nearly  all  the  energy  is  in 
the  air.  If  an  absorbent  strip  is  stuck  to  the  ferrite,  there  will 
be  a  large  attenuation  in  the  direction  of  propagation  corresponding 
to  (b)  and  a  low  attenuation  in  the  reverse  direction. 

It  is  very  difficult  to  calculate  field  displacement  sys¬ 
tems  because  of  the  complex  permeability  tensors.  A  relatively 
simple  case  is  that  of  a  thin  slab  in  the  plane  of  circular  polari¬ 
zation  of  the  guide:  in  this  case  the  permeability  would  be 
in  the  propagation  direction  and  p'  in  the  reverse  direction. 

We  can  have  either  p'"*^  <  p'  or  p'"*"  >  p'  depending  on 
the  value  of  H^.  It  is  essential  to  operate  in  the  former  range  so 
as  to  have  the  maximum  difference  between  the  two  directions  of 
propagation,  which  means  using  magnetic  fields  below  the  resonant 
field.  Under  these  conditions,  however,  the  system  is  only  of 


interest  at  high  frequencies  where  the  low-field  waves  are  negligi¬ 


ble. 
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2.4  The  Junction  Circulator 

A  particularly  compact  device  which  has  found  a  great 
many  applications  over  a  wide  frequency  range  Is  the  symmetrical 
Y-junctlon  circulator  as  shown  In  Fig.  (2.6). 


Fig.  (2.6)  (a)  The  circuit  symbol  for  the  Y-junctlon 
circulator  Indicates  transmission  as 
and  so  on, 

(b)  The  waveguide  version  show?  a  vertically 
magnetized  ferrite  rod  as  the  nonrecipro¬ 
cal  element, 

(c)  The  stripline  version  shows  the  ferrite 
disc  which  is  magnetized  perpendicular 
to  the  broad  face. 

The  feature  that  makes  this  device  so  attractive  is  its 
simple,  compact,  light  construction  at  all  frequencies.  Although 
it  has  been  designed  to  operate  in  bands  from  100  MHz  through  the 
ITHF  region  and  the  microwave  region  to  the  mm-wave  range,  it 


rarely  weighs  more  than  a  pound,  and  the  largest  models  are  less 
than  6  inches  in  diameter.  In  the  stripline  version  Fig.  (2.6(c)), 
it  is  often  manufactured  by  printed  circuit  techniques  and  is  minia¬ 
turized  through  the  use  of  high  dielectric  constant  materials,  which 
also  permits  it  to  carry  higher  r-f  power  before  corona  effects  set 
in.  The  device  will  tolerate  high  signal  power  considerably  in 
excess  of  50  Kw  (peak)  and  500  watts  (average)  when  operated  as  a 
circulator,  switch,  or  modulator. 

Until  very  recently,  the  Y  circulator  and  switches  have 
been  developed  empirically  and  explained  on  a  phenomenological 
basis  in  terms  of  the  scattering  matrix  treatment. 

The  circulation  mechanism  can  be  explained  briefly,  as 
follows:  The  wave  incident  at  Port  I  of  Fig.  (2.7)  is  the  fun¬ 
damental  mode  with  the  r-f  electric  vector  polarized  parallel  to 
the  direction  of  the  applied  magnetic  field.  Thus,  the  ferrite 
disc  supports  two  typical  electromagnetic  modes  characteristic  of 
such  a  tensor  magnetic  medium.  Two  mutually  orthogonal  components 
of  the  r-f  magnetic  field  are  excited  in  the  plane  perpendicular 
to  the  d-c  magnetic  field  vector,  i.e.,  the  plane  parallel  to  the 
broad  face  of  the  disc.  One  mode  has  a  right-hand  circularly 
polarized  r-f  magnetic  vector  and  has  increasing  field  intensities 
toward  the  right  side  of  the  input  port.  The  other  mode  has  the 
opposite  sense  of  rotation.  The  combination  of  these  two  modes 
has  a  hyperbolic  transverse  electric  field  dependence.  The  plane 
wave  incident  at  the  input  stripline  port  is  distorted  within  the 
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ferrite  so  that  field  Intensities  increase  toward  the  region  between 
input  Port  I  and  output  Port  II,  as  shown  by  the  contours  of  con¬ 
stant  electric  field  intensity  Fig.  (2.7). 


(Input) 


Fig.  (2.7)  Idealized  sketch  of  the  contours  of 

constant  r-f  electric  field  intensity 
within  the  ferrite  disc  of  the  stripline 
Y-circulator  illustrating  the  nonrecipro¬ 
cal  field  displacement  principle  of  circu¬ 
lation. 


This  combination  of  the  two  sets  of  rotating  modes  com¬ 
puted  by  Bosma^^^^  produces  a  wave  pattern  which  has  its  field- 
displacement  peak  between  the  input  and  output  ports  and  a  null  at 
the  Isolated  port.  Both  the  electric  and  magnetic  field  configura¬ 
tions  established  by  the  superposition  of  these  two  modes  result 
in  a  pattern  for  which  Poynting  vectors  exist  at  input  Port  I  and 
output  Port  II  and  are  oppositely  directed,  whereas  the  vector 


vanishes  at  Port  III. 
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2.5  Comparison  of  Isolator  Types 

In  this  chapter  four  basic  types  of  isolators  have  been 
discussed  with  respect  to  method  of  operation,  construction,  and 
design  consideration.  It  is  appropriate  to  compare  the  types  to 
bring  out  the  advantages  of  each.  Fig.  (2.8)  is  convenient  for 
this  purpose. 

The  field  displacement  isolator  requires  a  medium  biasing 
field  and,  thus,  a  medium  weight  magnet.  The  Faraday  rotation 
isolator  requires  the  lowest  biasing  field  and  seems  capable  of 
providing  isolation  approaching  those  of  the  field  displacement 
isolator.  Its  bandwidth  capabilities  are  somewhat  greater  than  the 
field  displacement  device  but  not  as  great  as  the  resonance  isola¬ 
tor.  The  power  handling  capabilities  of  both  the  field  displace¬ 
ment  and  Faraday  rotation  isolators  are  limited  by  tlie  loss  films 
employed  which  can  dissipate  a  relatively  small  amount  of  power. 

It  should  be  pointed  out  that  a  circulator  can  be  used 
as  an  isolator  by  terminating  the  unused  ports.  This  can  be  par¬ 
ticularly  advantageous  from  a  power  handling  standpoint,  since 
reflected  power  can  be  dissipated  in  a  termination  which  is  capable 
of  being  cooled,  rather  than  in  the  f errimagnetic  material.  The 
Y-junction  circulator,  when  used  as  an  Isolator,  rates  an  E  on  all 
counts  in  Fig.  (2.8). 


Fig.  (2.8)  Comparison  of  isolator  types . 


CHAPTER  3 


ANALYTICAL  FORMULATION 

This  chapter  deals  mainly  with  explaininj.'  tie-  derivation 
of  tile  analytical  expression  of  tlie  proposed  structure.  I'lie  struc¬ 
ture  (Fig.  3-1)  consists  of  a  dielectric  waveguide  containing  a 
ferrite  layer  in  which  periodic  metal  strips  are  provided.  This 
new  structure  is  an  Improvement  of  the  recently  proposed  device  in 
which  a  notched  grating  was  used  Instead  of  the  metal  strips^ 

In  thi’  analysis  it  has  been  assumed  that  the  structure  is  infinitely 
periodic  in  the  x-direction  and  that  it  is  invariant  in  the  y-direc- 
tions.  The  metal  strips  are  also  assumed  to  be  infinitely  thin. 

The  operating  principle  of  this  device  lies  in  the  fact 
that  the  propagation  constants  of  the  guided  waves  in  the  opposite 
directions  are  different  when  a  dc  magnetic  field  is  applied.  When 
a  periodic  perturbation  is  provided  along  the  waveguide,  it  is  pos¬ 
sible  at  a  certain  frequency  interval  that  a  guided  wave  propagating 
in  one  direction  becomes  leaky  and,  hence,  it  attenuates  as  it 
propagates;  whereas,  the  one  propagating  in  the  opposite  direction 
is  still  guided  without  leakage.  Hence,  the  structure  works  as  an 
i  solator . 

In  the  structure  proposed  previously,  the  periodic  per¬ 
turbation  was  produced  by  a  grooved  grating.  It  turns  out  that  the 
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grooved  grating  provides  only  a  small  perturbation*  .  In  the  new 
structure  the  nu-tal  strips  would  provide  an  appreciable  amount  of 
perturbation,  since  the  principal  electric  field  is  chosen  to  be 
parallel  to  the  strips.  Hence,  the  attenuation  constant  due  to 
leakage  in  this  new  system  would  also  be  large  and  the  entire  device 
length  can  be  reduced  for  a  given  degree  of  isolation. 

In  the  analysis  it  is  assumed  that  the  problem  in  invari¬ 
ant  in  the  y-direction.  It  is  also  assumed  that  the  strips  have 
no  thickness  and  only  E^,  and  components  are  present.  The 
problem  is  first  solved  for  the  system  without  metal  strips  using 
the  transverse  resonance  technique.  Then  by  applying  Floquet's 
theorem,  the  electric  fields  outside  the  strip  and  the  currents 
in  the  strips  are  represented  in  terms  of  space  harmonics.  The 
transverse  equivalent  circuit  of  each  space  harmonic  is  then 
derived  with  the  effect  of  the  metal  strip.  The  next  step  is  the 
application  of  Galerkin's  method  to  the  entire  spectrum.  In  this 
method  the  current  is  expanded  in  terms  of  known  basis  functions. 

The  resulting  equation  could  then  be  used  to  solve  for  the  propa¬ 
gation  constant  of  the  grating  structure. 


Metal  Strips 


3.1 


52 


Tht'  Transfer  Matrices  and  the  Equivalent  Transverse 
Resonance  Circuit 

f  19  ] 

The  transverse  resonance  method  is  a  powerful  way  to 

derive  a  transcendental  equation  in  terms  of  the  propagation  con¬ 
stant  of  the  system  in  question.  The  sum  of  the  impedances  at  one 
point  in  the  transverse  direction  is  zero  (resonance)  only  for  the 
correct  values  of  propagation  constant. 

In  the  case  of  a  multi-layered  structure,  formulation  of 
the  impedance  in  the  transverse  direction  can  be  facilitated  by  the 
use  of  transfer  matrices  for  each  layer  in  question.  In  this  way 
one  could  easily  obtain  the  impedance  at  any  plane  in  the  system 
with  respect  to  the  impedance  at  a  given  reference  plane. 

The  derivation  of  the  transfer  matrix  for  a  ferrite  layer 
is  as  follows.  Consider  the  structure  in  Fig.  3.2. 


Fig.  (3.2). 


From  Maxwell's  equations 


VxE  =  -j  (Op  H 


(3 


For  a  TF.  wave  in  the  x  direction  we  are  going  to  have  E  ,  11  and 

°  y  X 

components  only. 
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Equatin  for  a  and  a  in  Maxwell's  equations  and  solving  for  H  ( 

X  Z  X 

we  get 

H  (z)  =  — ^ -  ta(a6  +  j  n)  e  +  b(a6  -  j  n)  (3 


X  up  P 

o  e 


2  2  2  ? 
where  p  =(p  -  K')/p  ,  a=K/p  ,u=l+tuu/(u  -u)] 

e  o  n  o 
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2  2  /  2  ^ 
K  =  WO)  /  (ai  -  w  )  ,  u  =  yH  ,  ui  =  y47rM  and  n=/w  yijce  -  B 

no  o  m  .jooeom 

Here,  y  Is  the  gyromagnetic  ratio  and  4ttM  is  the  saturation  magneti¬ 
zation.  Also,  tile  time  and  propagation  factor  exp  { j  (ait  -  Bx)  }  is 
to  be  understood  and  suppressed. 

We  will  now  Identify  E  and  H  as  V  and  —I,  respectively. 
Then  the  values  of  V (o)  and  I(o)  at  Z  =  0  are  related  to  those  at 
V(z)  and  I(z)  at  z  =  z  by 


cos  nz  - 


sin  r\z 


wp  u 
o  e 

j - sin  nz 

n 


b  nz  +  ^  sin  nzj  [l(z)J 

°  ®  (3.7) 

This  cxperssion  relates  the  impedance  at  one  point  in  the  transverse 

direction  to  a  reference  value  of  the  impedance. 

In  the  case  of  a  dl'^lectric  material,  the  transfer  matrix 

can  be  written  as 


cos  nz 


i  -  sin  nz  V(z) 

n 


1  -  sin  nz  cos  nz 

wp 

.  o 


(3.8) 


This  can  be  derived  by  letting  o  =  0  and  Pp=  1  in  (3.7) 

The  transverse  resonance  equivalent  network  for  the  structure  in 

Fig.  3.1  is  given  by  removing  the  current  source  in  Fig.  3.3.  The 
current  source  appears  when  the  equivalent  transverse  network  con¬ 
cept  is  Introduced  in  the  analysis  of  the  grating  structure  with 
metal  strips.  Such  a  process  will  be  explained  in  the  next  section. 


I 
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The  eigenvalue  equation  +  Yj^  =  0  provides  the  propagation  con¬ 
stant  B.  Note  that  we  obtain  two  different  values  of  B  corresponding 
to  waves  propagating  in  the  opposite  directions. 


z  = 


Z  =0 


z  =  -t. 


Air-  Region 


Ferrite-  Slob 


t  ^  Jn 

I  YL(^n)  "'"Yu 


Dielectric  -  Slob 


Air  -  Region 


Fig.  (3.3)  Transverse  equivalent  network 
for  the  space  harmonics. 
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3.2  Floquet's  Theorem  and  the  Space  Harmonics 


The  basis  for  a  study  of  periodic  transmission  systems  Is 
a  theorem  due  to  a  French  mathematician,  Floquet.  This  theorem  may 
he  stated  as  follows:  For  a  given  mode  of  propagation  at  a  given 
steady-state  frequency,  the  fields  at  one  cross  section  differ  from 
those  a  period  away  only  by  a  complex  constant.  The  theorem  is 
true  whether  or  not  the  structure  contains  loss  as  long  as  it  is 
periodic.  The  proof  of  the  theorem  lies  in  the  fact  that  when  the 
structure  having  infinite  length  is  displaced  along  its  axis  by  one 
period,  it  cannot  be  distinguished  from  its  original  structure. 
Hence,  we  write  the  electric  or  magnetic  field  as 


f.  .  R(,,  y,  .) 


(3.9) 


where  E(x,  y,  z)  is  periodic  in  x  with  period  d. 

—  Yx 

The  field  F.(x,  y,  z)  e  can  be  expanded  in  a  Fourier 
series  of  the  form 


K(x,  y,  7.)  e 


-YX 


To  find  F.^(y,z),  multiply  the  above  by  e 


e-"’^  (3.10) 

all  n 

j  (2'mii/d)x  •+  YX 


E(x,  y,  z)  E^(y,  zl  e 


j(2'nm/d)x  -  j(2Tm/d)x 


all  n 


(3.11) 
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Integrate  both  sides  from  to  +  d 


.Xj+d 

I  f'-  (x ,  y ,  Z 

J  V 


.  j (?im/d)x 
^)  c-^  dx  = 


eT 

all  n  J  X, 


E^(y,z)  e 


j  (27/d)  (ni-n)x 


dx 


=  0  ,  in  ^  n 

=  E^(y,z)d  ,  m  =  n 


(3.12) 


Therefore, 


E 

n 


(y,z) 


1 

d 


[E(x,y,z)  e 


[y+j  (2Tm/d)]x  , 
e  dx 


(3.13) 


This  expression  is  used  later  to  calculate  the  amplitude  of  indivi¬ 
dual  terms  in  the  series. 

The  n-th  term  on  the  right  of  the  series  representation 
is  called  the  n-th  space  harmonic.  It  has  a  propagation  constant 
Y  +  j  If  there  are  no  losses  in  the  system,  it  turns  out  that 

y  is  either  purely  real  or  purely  imaginary.  If  imaginary,  we 
write  Y  =  j  p.^  and  define 

P  =  6  +  ^  (3.1A) 

nod 


where  6^  is  called  the  phase  constant  for  the  n-th  space  harmonic. 
A  velocity  v^  -  phase  velocity  of  the  n-th  space  har¬ 
monic  . 
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Applying  Floqiiet's  theorem  to  the  structure  containing 
periodic  metal  strips,  we  can  write  expressions  for  the  electric 
field  and  the  current  on  the  strips  as 

jB  X 

•e  "  (3.15) 


CO 


J 

y 


j 

n 


JS„x 


6  -  6  +  -^  (3.16) 

n  d 


Since  and  are  expressed  in  the  same  manner,  we  only 

need  to  relate  E  and  J  .  Since  E  and  J  are  no  longer  functions 
n  n  n  n  ^ 

of  X,  we  can  use  the  spectral  domain  immittance  method.  For  each 
space  harmonic,  each  layer  is  replaced  with  a  transfer  matrix. 

From  the  transverse  equivalent  circuit  (Fig.  3.3)  the  space  har¬ 
monics  of  the  E  at  z  =  0  and  of  J  are  related  by 

y 


^n  ~  Y„(6  )  +  Y,  (6  ) 


U'  n' 


L  n 


J  =  Z  J 
n  n  n 


(3.17) 


The  quantities  Y  and  Y  are  input  admittances  looking  upward  and 

U  L* 

downward  at  z  =  0.  They  could  easily  be  found  by  the  use  of  the 
transfer  matrix  methods  discussed  in  the  previous  section.  In 
fact,  it  is  only  necessary  to  substitute  8^  for  3  in  the  expres¬ 


sions  derived  in  that  section. 
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This  process  is  identified  as  Galerkin's  procedure.  The  above 
equation  has  a  unique  solution  only  if 

det  [K..]  =  0  (3.22) 

ij 

From  this  equation,  the  dispersion  characteristics  can  be  obtained, 
and  we  can  compute  6  for  a  given  frequency. 


CHAPTER  4 


THE  COMPUTER  PROGRAM 

A  rompdter  program  has  been  written  to  analyze  the 
ferritt-diclectric  structure  containing  metal  strips  (Fig.  3.1). 
The  program  is  based  on  the  theory  discussed  in  Chapter  3.  The 
best  way  to  explain  how  the  theory  has  been  implemented  numeri¬ 
cally  by  the  program  is  to  look  at  the  flow  chart  of  the  program 
(Fig.  4.1)  relating  how  the  subroutines  are  linked.  A  print-out 
of  the  program  is  given  in  the  Appendix. 


! 

i 
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Fig.  (4.1). 


— .t' 
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.  1  Program  PERLOAR 

This  is  the  main  program  which  controls  all  of  the  sub¬ 
routines.  First  the  data  defining  the  structure  dimensions  and 
material  properties  are  read  and  then  printed.  The  space  har¬ 
monics  are  then  calculated  and  stored  in  an  array  by  the  subroutine 
SPACE. 

Polders  permeability  tensor  is  then  calculated,  Oreen^^®) 
and  Sandy's  formulas  are  used  for  the  ferrite  below  saturation. 

The  program  then  calculates  the  propagation  constant  of 
the  structure  without  periodic  metal  strip.t  using  Subroutine 
UNIFORM.  This  result  is  then  used  as  a  first  approximation  of  B 
for  the  structure  with  the  metal  strips  where  the  actual  theory 
discussed  in  Chapter  3  is  implemented. 

Subroutine  PERIOD  is  then  called.  This  suliroutine  in 
turn  calls  the  rest  of  the  subroutines  to  implement  tlie  root  search 
of  the  actual  structure.  This  will  be  discussed  in  more  detail  in 
the  next  section.  This  subroutine  will  give  the  result  of  3  for 
the  re.spective  frequencies  desired.  The  program  then  calls  sub¬ 
routine  GRAPH;  this  is  a  ZETA  plot  routine  to  plot  the  K-3  diagram 


desired . 
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4.3  Subroutine  UNIFORM 

This  routine  searches  for  a  solution  of  6  when  the 
structure  consists  of  just  a  dielectric-ferrite  guide.  Subroutine 
UNIFORM  is  called  to  give  the  impedance  looking  up  and  down  trans¬ 
versely  at  the  ferrite  dielectric  boundary.  According  to  trans¬ 
verse  resonance  the  solution  of  6  occurs  when  the  sum  of  the  jjnped- 
ance  looking  up  and  looking  down  is  zero. 

This  routine  searches  only  for  the  real  part  of  B.  The 
way  this  is  done  is  by  fixing  a  lower  and  upper  bound  where  6  may 
exist  and  then  scan  this  region  by  dividing  it  into  100  segments. 
The  lower  bound  has  been  chosen  to  be  the  6  of  air  only  and  the 
upper  bound  the  6  for  ferrite  only,  since  ferrite  usually  has  the 
largest  value  of  6.  The  only  problem  that  may  occur  in  this  rou¬ 
tine  is  when  the  effective  permeability  of  the  ferrite  becomes 
negative.  In  this  case  either  the  value  of  B  could  be  smaller  than 
that  of  air  or  no  solution  would  be  found  as  has  been  frequently 
experienced . 

4.4  Subroutine  PERIOD 

This  routine  is  mainly  a  root  search  routine.  The  func¬ 
tion  is  mainly  to  find  B  which  makes  the  determinant  zero.  The 
determinant  is  calculated  by  subroutine  DETER,  so  this  routine 
gives  the  value  of  B,  plugs  it  in  DETER  and  then  tries  to  find  the 
value  of  B  for  zero  determinant. 
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Muller's  method  is  used  to  search  for  3.  Muller's  method 
is  an  improvement  on  Newton's  method.  For  more  detail  one  should 
refer  to  Peter  Henrici,  "Elements  of  Numerical  Analysis,"  John 
Wiley  &  Sons,  1964,  page  198. 

This  routine  could  be  replaced  hy  an  IMSL  routin  such  as 
ZANALY.  ZANALY  has  been  used  and  tested;  it  actually  gives  the 
same  results  as  subroutine  PERIOD.  The  only  difference  is  that 
using  IMSL  is  more  expensive  and  one  has  less  control  of  what  goes 
on. 


4,5  Subroutine  DETER 


The  function  of  this  routine  is  to  calculate  the  deter¬ 
minant  for  a  given  input  value  of  3.  First,  the  matrix  is 
formed  by  putting  together  the  space  harmonics  of  the  basis  func¬ 
tions.  It  has  been  found  experimentally  that  taking  ±40  space  har¬ 
monics  would  give  an  accurate  answer.  In  order  to  calculate  the 
symmetric  matrix  the  routine  calls  subroutine  IMPED  which  in 

turn  calls  subroutine  EQUIV  to  get  the  necessary  information  on  the 
impedance  looking  in  the  upper  and  lower  regions,  to  be  able  to 
implement  the  formula 


K.  .  = 
ij 


rY,,(3  )  +  Y  (3  )] 
n  |_  U  n  L  n  J 


(4.4) 


The  determinant  is  then  calculated  using  standard  equations. 


CHAPTER  5 


THE  EXPERIMENTAL  INVESTIGATION 

This  chapter  deals  mainly  with  the  experiments  that  have 
been  carried  out.  Both  theoretical  and  experimental  results  are 
included.  The  main  purpose  of  this  investigation  is  to  try  to  see 
if  the  numerical  results  are  reasonable.  Another  purpose  is  to 
investigate  experimentally  if  this  structure  could  be  used  as  an 
isolator . 

Extensive  experiments  have  been  performed  to  investigate 
the  validity  of  the  theory.  Since  the  theory  is  based  on  the 
assumption  that  the  structure  is  infinite  in  both  X  and  Y  direc¬ 
tions,  and  practically  no  such  device  could  be  constructed,  an 
experimental  investigation  has  been  carried  out  to  sec  the  effect 
of  finiteness  of  the  structure. 
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5.1  Evaluation  of  the  Structure  for  Experimental  Purposes 

The  structure  analyzed  in  theory  could  he  simulated  by  two 
kinds  of  guides,  n£imely  dielectric  guides  or  image  guides.  Both 
systems  have  their  advantages  and  disadvantages,  but  the  fact  that 
it  is  difficult  to  support  dielectric  guides  makes  the  image  guide 
more  feasible  to  build  an  Integrated  system.  Although  experiments 
have  been  carried  out  using  both  systems,  the  image  guide  system  has 
been  finally  adopted  for  these  experiments.  A  brief  description  of 
the  two  systems  used  is  given  below. 

(a)  Dielectric  Guide  System 

The  structure  used  was  made  of  Stycast  Hi-K  (e^  =  10)  as 
the  dielectric  and  of  Magnesium  ferrite.  The  dielectric  guide  was 
fed  from  two  identical  horns.  The  horns  have  waveguide-to-coaxial 
adapters  and  the  system  was  tested  using  a  coaxial  ref lactometer . 

A  schematic  diagram  of  the  entire  system  is  given  in  Fig.  5.1.  The 
setup  was  mechanically  not  very  stable  because  the  positioning  of 
the  dielectric  between  the  horns  was  very  sensitive  to  small  move¬ 
ments,  Another  major  problem  with  this  system  was  that  the  reflec¬ 
tion  coefficient  could  not  be  well  monitored  since  the  coaxial-to- 
waveguide  adapter  produced  significant  reflections. 

Cb)  Image  Guide  System 

Two  systems  using  image  guides  have  been  considered.  One 
system  used  the  coaxial  ref lactometer  just  like  the  dielectric  guide 


system,  and  the  other  system  was  designed  to  work  with  a  waveguide 
ref lectometer.  The  use  of  the  waveguide  ref lactometer  was  intended 
to  eliminate  the  reflections  caused  by  the  waveguide  to  coaxial 
transition.  A  diagram  explaining  the  image  guide  system  is  shown 
in  Fig.  5.2.  The  waveguide  ref lectometer  had  a  return  loss  of  about 
-40  dB  while  the  coaxial  ref lectometer  with  horn  to  horn  position 
had  a  return  loss  of  -15  dB. 

A  problem  with  the  image  guide  system  as  compared  to  the 
dielectric  guide  system  is  whether  or  not  the  metal  strips  should 
contact  the  ground  plane.  This  problem  will  be  investigated  in 
detail  in  a  later  section. 

5.2  The  Test  Setup 

The  final  choice  of  the  setup  was  the  image  guide  system 
using  the  waveguide  ref lectometer .  Due  to  the  X-band  waveguides 
used,  the  frequency  swept  was  from  7  GHz  to  16  GHz.  The  measure¬ 
ments  taken  were  mainly  insertion  loss  and  return  loss  using  the 
H.P.  scalar  network  analyzer.  The  magnet  used  was  a  water  cooled 
Harvy  Wells  magnet  which  could  go  up  to  8  KG  with  the  supply  avail¬ 
able.  A  Gauss  meter  by  F.  W.  Bell  was  used  to  monitor  the  magnetic 
field.  A  diagram  of  the  whole  setup  is  given  in  Fig.  5.3. 


5.3  Preliininary  Investigation,  Non-reciprocal  Effects  on 

Applying  an  External  Magnetic  Field 

This  investigation  has  been  made  to  answer  some  of  the 
questions  that  could  be  put  forth  on  the  structure  used.  The  dielec¬ 
tric  used  was  Stycast  HiK  (e  =  10)  and  the  ferrite  was  Magnesium 
type  (TT-1-105) .  The  ferrite  and  dielectric  widths  were  3  mm  each 
and  the  height  of  the  image  guide  was  5  mm.  The  period  and  strip 
width  of  the  grating  were  13  and  3  mm,  respectively. 

(a)  No  Metal  Strips 

This  experiment  has  been  carried  out  to  verify  that  the 
waves  behave  in  a  reciprocal  manner  when  the  metal  strips  are  not 
inserted.  Another  purpose  of  this  experiment  is  to  show  at  what 
frequency  the  ferrite  resonates  for  a  particular  magnetization. 

From  the  results  of  the  experiments  shown  in  Fig.  5.4,  one  could 
easily  see  that  in  one  direction  the  system  goes  in  resonance  at  a 
lower  frequency  than  the  other  direction.  This  behavior  can  not 
be  explained  as  yet,  but  from  what  has  been  observed  the  difference 
in  frequency  is  only  200  MHz  at  most.  This  system  could  be  used 
for  isolation  purposes,  but  the  bandwidth  is  very  narrow. 

(b)  Metal  Strips  Touching  Ground  Plane 

This  experiment  has  been  performed  as  part  of  the  investi¬ 
gation  to  determine  whether  or  not  the  metal  strips  should  touch 
the  ground  plane.  From  the  experimental  results  of  Fig.  5.5,  one 
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Fig,  (5.Aa)  Guide  with  no  metal  strips. 
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Fig. (5,5b)  Metal  strips  touching  ground  plane. 
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could  see  that  for  no  magnetization  almost  all  the  energy  is 
scattered  or  attenuated.  As  the  magnetic  field  is  applied,  large 
non-reciprocity  is  evident  at  two  frequency  ranges,  namely  from  8 
to  9  GHz  and  11.5  to  12.5  GHz.  One  notes  that  this  occurs  at  a 
magnetization  of  only  about  .5  KG,  indicating  that  the  system  is 
way  below  saturation  and,  of  course,  resonance. 

This  structure  has  not  been  considered  further  since  at 
zero  magnetization  the  theory  does  not  predict  that  all  the  energy 
should  be  scattered.  Although  the  insertion  loss  in  the  forward 
direction  is  about  10  dB,  this  structure  may  be  useful  if  attenua¬ 
tion  is  reduced  by  optimization. 

(c)  Use  of  Low  Dielectric  Constant  Slabs  to  Suspend  the  Guide 

Above  the  Ground  Plane  so  that  the  Strips  would  not  Touch 

This  experiment  has  also  been  part  of  the  Investigation 
to  find  out  whether  the  strips  should  touch  the  ground  plane  or 
not.  In  the  experimental  results  shown  in  Figs.  5.6  and  5.7,  thin 
dielectric  films  of  e  =  2.3  and  0.7  mm  thickness  have  been  inserted 
between  the  ferrite-dielectric  rod  and  the  ground  plane.  From  the 
result,  it  is  evident  that  a  stopband  and  a  leaky  wave  region  could 
be  seen  at  8  and  9.5  GHz,  respectively. 

As  compared  to  when  the  strips  were  touching  the  ground 
plane,  the  non-reciprocity  at  12  GHz  is  still  evident  at  low  mag¬ 
netization  in  the  present  case.  The  other  phenomenon  is  that  in 
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Fig.  (5.7b)  Guide  suspended  1.4  mm  above  ground 
plane  by  means  of  dielectric 
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one  direction  the  system  goes  in  the  leaky  wave  region  at  a  different 
frequency  than  the  other  direction.  This  difference  is  only  100  MHz 
and  could  not  be  of  much  use.  One  also  notices  that  in  the  negative 
permeability  region  non-reciprocity  occurs.  The  bandwidth  of  the 
non-reciprocal  effect  is  about  1.5  GHz.  This  effect  is  very  visible 
especially  near  3  KG, 

To  find  out  the  effect  of  the  thickness  of  the  low  dielec¬ 
tric  constant  layer,  experiments  have  been  performed  for  two  differ¬ 
ent  values  of  thickness.  Fig.  5.6  is  for  one  slab  of  .7  mm  thickness 
while  the  data  of  Fig.  5.7  is  for  two  slabs  of  . 7  mm  thickness  each, 

A  comparison  of  the  results  of  Fig.  5.6  to  those  of  5.7,  shows  that 
the  overall  attenuation  decreases  with  increase  of  the  thickness  of 
the  low  dielectric  constant  layer.  On  the  other  hand,  the  non¬ 
reciprocity  and  the  stopband  are  larger  for  a  thinner  layer  of  low 
dielectric  constant  material. 

(d)  Use  of  Metal  Strips  which  do  not  Touch  Ground  Plane: 

Strips  are  truncated  at  1  mm  Above  Ground  Plane 

Another  way  of  preventing  the  metal  strips  from  touching 
the  ground  plane  is  to  truncate  the  metal  strips  on  the  ground  plane 
side.  The  results  obtained  with  this  structure  are  similar  to  the 
ones  when  the  low  dielectric  slab  was  used.  The  only  difference  is 
that  the  overall  attenuation  is  2  dB  less.  This  is  mainly  due  to 
the  use  of  shorter  metal  strips.  The  results  of  this  investigation 
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are  shown  In  Fig.  5.8.  The  stophand  in  this  system  is  hardly  visible 
at  8  GHz  since  the  metal  perturbation  is  small, 

(e)  The  Effect  of  the  Period  and  Strip  Width 

The  change  in  period  of  the  metal  strips  should  only  change 
the  position  of  the  stopband,  although  it  may  have  an  effect  on  the 
bandwidth  of  the  non-reciprocal  part  of  the  result.  Experimentally 
not  much  has  been  made  in  this  respect.  The  only  experiment  per¬ 
formed  is  with  the  period  and  strip  width  11  mm  and  1.5  ram,  respec¬ 
tively  instead  of  13  mm  and  3  mm  as  in  the  previous  case.  The 
results  of  this  experiment  are  shown  in  Fig.  5.9. 

From  the  results,  one  could  see  that  the  overall  attenua¬ 
tion  is  reduced  compared  to  Fig.  5.8,  while  the  non-reciprocal 
effects  are  still  preserved.  This  observation  makes  sense  since  the 
narrower  the  metal  strips,  the  less  scatter  and  perturbation  occurs. 

5.4  Effect  of  having  a  Finite  Grating  Structure 

In  this  sect'' on,  we  study  the  effect  of  finite  length  of 
the  grating  section.  The  analysis  of  a  periodic  structure  applies 
to  an  infinitely  long  structure  in  the  strict  sense.  It  is  expected 
that  the  characteristics  of  a  grating  of  finite  length  will  differ 
somewhat  from  those  of  an  infinitely  long  one.  We  will  study  such 
difference. 

The  structure  used  for  this  Investigation  was  made  of 
Alumina  (e  %  10),  and  no  ferrite  has  been  used.  The  structure  was 
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made  with  alumina  on  both  sides  of  the  metal  strips.  Two  alumina 
slabs  were  glued  together  with  the  metal  strips  In  the  middle.  The 
period  and  strip  width  was  IQ  mm  and  2  mm,  respectively.  The  thick¬ 
ness  of  the  alumina  slabs  was  1.5  mm  and  the  height  4  mm. 

The  experiment  consisted  of  starting  with  17  metal  strips 
and  then  taking  off  four  strips  at  a  time  to  measure  the  system  with 
13,  9,  5,  3,  1  and  finally  no  metal  strips  at  all.  In  this  experi¬ 
ment,  the  metal  strips  were  isolated  from  the  ground  plane.  Each 
of  the  seven  structures  described  above  were  tested  with  five  dif¬ 
ferent  types  of  Isolation  from  the  ground  plane.  In  one  of  these 
five  methods,  the  structure  was  isolated  from  the  ground  plane  with 
a  layer  of  3M  scotch  tape.  In  the  remaining  four  cases,  four  slabs 
of  dielectric  material  (e  =  2.3)  of  .7  mm  thickness  each  were 
placed  under  the  alumina  structure  one  by  one  to  make  different 
thicknesses  of  low  dielectric  constant  material  under  the  main 
guide.  The  experimental  results  of  this  investigation  are  shown 
in  Fig.  5.10. 

The  computer  program  has  been  run  for  these  particular 
cases  to  investigate  if  there  is  any  correlation  with  the  experi¬ 
mental  results.  The  computer  result  is  shown  in  Fig.  5.11,  where 
the  attenuation  constant  is  given  in  Nepers  per  mm  x  (-j)  . 

The  computer  results  show  that  the  system  should  have  a 
stopband  between  10.5  GHz  and  13  GHz.  From  the  experimental  result, 
the  stopband  is  between  10  and  11.5  GHz,  so  there  is  a  discrepancy 
of  about  1  GHz  in  the  bandwidth  of  the  stopband.  The  overall  shape 


metal  strips.  This  experiment  shows  the  effect  of  varyii 
the  thickness  of  a  low  dielectric  constant  material  (C“ 
placed  between  the  alumina  and  the  ground  plane  to 
Isolate  the  metal  strips. 
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of  the  theoretical  and  exper Jjsental  results  are  very  similar.  It 
is  evident  theoretically  that  the  attenuation  at  frequencies  below 
the  stopband  Is  negligible.  Experimentally,  the  insertion  loss 
curves  for  the  cases  with  and  without  metal  strips  differ  very  lit¬ 
tle  for  the  region  below  the  stopband.  This  confirms  negligible 
loss  predicted  by  theory.  In  the  leaky  region  (which  occurs  at 
14  GHz  theoretically  and  at  12  GHz  experimentally) ,  the  attenua¬ 
tion  predicted  is  about  1  dB/cm.  This  ratio  gives  us  about  10  dB 
of  Insertion  loss  for  the  particular  length  of  device  used  in  the 
experiment.  The  experimentally  observed  insertion  loss  is  also 

about  10  dB.  For  this  particular  case,  the  computer  program  has 

125] 

been  run  using  the  effective  dielectric  constant  method  to 
compensate  for  making  the  structure  finite  in  the  y  direction. 

The  effective  dielectric  constant  has  been  obtained  in 
the  following  manner.  The  propagation  constant  of  a  dielectric 
guide  with  a  height  of  100  mm  was  first  computed.  We  now  assume 
that  this  guide  has  an  infinite  height  ana,  hence,  is  a  slab  wave¬ 
guide.  Next,  we  adjust  the  value  of  the  dielectric  constant  of  a 
guide  with  a  height  of  4  mm  in  such  a  way  that  the  propagation 
constant  is  identical  to  the  one  in  the  guide  with  an  "Infinite" 
height.  Such  a  value  of  dielectric  constant  is  taken  as  the 
effective  dielectric  constant  of  the  4  mm  guide.  Some  standard 
software  to  calculate  the  propagation  constant  of  an  image  guide 
has  been  used  to  calculate  the  effective  dielectric  constant.  The 
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program  has  been  run  for  each  particular  frequency  since  the 
effective  dielectric  constant  is  very  dispersive,  especially  for 
a  small  guide  like  the  one  used  here. 

One  should  note  that  the  effective  dielectric  constant 
method  only  takes  care  of  a  finite  size  dielectric,  but  does  not 
take  care  of  finite  length  of  metal  strips.  On  the  other  hand,  it 
is  very  difficult  to  do  the  same  effective  dielectric  constant 
method  for  the  case  where  ferrite  is  involved  since  no  extensive 
study  is  available  to  date  to  calculate  the  effective  dielectric 
constant  and  the  effective  permeability  of  a  ferrite  guide 

5.5  Effect  of  Making  the  Structure  Finite  in  the  Y-Direction 

It  was  pointed  out  in  the  previous  section  that  the 
effective  dielectric  constant  method  cannot  be  applied  to  the 
structure  containing  ferrite.  It  is,  therefore,  necessary  to 
experimentally  study  the  effect  of  a  finite  size  structure.  It  is 
inferred  that  the  waveguide  with  a  larger  dimension  exhibits  better 
agreement  with  the  theory  based  on  the  infinitely  high  structure. 

To  make  this  investigation,  three  different  sizes  of 
guides  were  produced.  Each  of  the  three  guides  had  a  different 
height.  One  was  5  mm  in  height  and  the  others  were  in  turn  10  mm 
and  15  mm.  The  period  of  the  metal  strips  used  was  10  ram  and  the 
strip  width  was  2  ram.  On  the  basis  of  the  previous  investigation 
(.see  5.4),  the  number  of  strips  used  was  11.  The  thickness  of 
ferrite  was  1.5  mm  and  the  thickness  of  dielectric  also  1.5  mm. 
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The  dielectric  was  Stycast  HIK  (e^  =  IQ)  and  the  ferrite  material 
had  a  dielectric  constant  of  12  and  a  saturation  magnetization  of 
1.75  KG.  Each  of  the  three  structures  was  run  in  the  magnet  with 
and  without  metal  strips.  One  run  was  made  with  just  one  layer  of 
3M  scotch  tape  and  one  with  one  slab  of  dielectric  of  . 7  mm  thick¬ 
ness  and  2.3  dielectric  constant.  Experiments  have  also  been  made 
with  a  dielectric  material  only  on  both  sides  of  the  metal  strips. 

The  experimental  results  of  this  Investigation  are  shown 
in  Fig.  5.12.  Figs.  5.13  and  5.14  are  tables  listing  the  observa¬ 
tions  worth  noticing  from  these  experimental  results.  The  most 
Important  observation  from  those  results  is  that  non-reciprocity 
occurs  in  the  negative  permeability  region.  The  best  example  of 
this  is  shown  at  a  magnetization  of  3  K  Gauss.  In  the  experimental 
results,  the  non-reciprocity  is  seen  clearly  for  10  mm  and  15  mm 
thickness.  The  region  of  non-reciprocity  occurs  from  11  GHz  to 
12  GHz.  This  frequency  range  corresponds  directly  with  the  nega¬ 
tive  permeability  region.  Although  the  negative  permeability 
region  is  from  10.6  GHz  to  13.3  GHz,  the  permeability  is  mostly  a 
large  negative  number  between  11  GHz  and  12  GHz.  Fig.  5.15  shows 
how  theoretically  the  negative  permeability  region  occurs. 

The  comouter  program  has  been  used  to  investigate  if 
this  non-reciprocal  effect  could  be  predicted  theoretically. 
Unfortunately,  for  the  structure  with  no  metal  strips  in  the  posi¬ 
tive  direction,  no  solution  was  found  for  the  negative  permeability 
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f'ig»(5.14)  Comparison  of  the  results  when  using  various  guide 
heights,  for  the  structure  with  metal  strips. 
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region  between  11  and  12  GHz.  Therefore,  no  initial  guess  was 
available  to  determine  how  the  structure  behaves  with  metal  strips. 

From  the  theoretical  results,  it  is  evident  that  more  than 
one  mode  exists  in  the  system.  From  the  computer  results  at  3  KG 
(Fig.  5.16),  it  could  Be  easily  seen  that  two  modes  exist  for  the 
system  both  with  and  without  metal  strips.  These  two  modes  may  be 
identified  as  the  magnetostatic  mode  and  the  dynamic  mode. 

Most  probably  the  non— reciprocity  occurring  in  the  nega¬ 
tive  permeability  region  is  due  to  the  dynamic  mode,  since  according 

f7] 

to  the  work  done  by  Bolle,  this  mode  is  the  most  dominant  in  this 
region . 

So,  since  no  initial  guess  was  available,  varluos  arbitrary 
initial  guesses  were  run  through  the  computer  program  for  the  deter¬ 
mination  of  the  propagation  constant.  Fig.  5.17  and  Fig.  5.18  show 
the  results  obtained  in  the  frequency  range  from  10  to  13  GHz  with 
different  values  of  initial  guess.  It  is  evident  that  when  a  small 
initial  guess  was  chosen,  there  was  no  non-reciprocity  in  the  attenua¬ 
tion  coefficient.  On  the  other  hand,  when  a  large  initial  guess  was 
used,  non-reciprocity  is  evident  in  the  11  to  12  GHz  frequency 
range.  From  the  computer  result,  the  forward  loss  should  be  about 
30  dB  while  the  backward  loss  should  be  about  60  dB.  Experimentally, 
with  a  structure  of  10  ram  height,  the  forward  loss  is  about  15  dB 
while  the  backward  loss  is  about  38  dB.  Although  these  results  do 
not  correlate  exactly,  one  could  infer  what  takes  place  since  the 
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greater  the  height  of  the  structure,  the  closer  the  correlation  with 
the  theory. 
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CHAPTER  6 


CONCLUSIONS 

6.1  The  Theoretical  Formulation 

The  theory  used  in  this  article  assumed  that  the  structure 
is  infinite  in  the  x  and  y  directions.  Having  the  structure  finite 
in  the  x-direction  may  not  he  very  far  from  Being  infinite  since  all 
the  parameters  derived  from  the  program  are  per  unit  length  in  the  x 
direction.  When  the  number  of  strips  is  finite,  one  must  then  con¬ 
sider  the  edge  effects  due  to  the  first  and  last  strips.  For  an 
appreciable  length  of  the  periodic  structure,  the  edge  effect  may, 
however,  be  neglected. 

The  main  deficiency  of  the  theory  is  the  assumption  that 
the  structure  is  infinite  in  the  y-direction.  Development  of  the 
theory  for  a  finite  y-dlmension  is  very  difficult,  since  one  has 
also  to  take  into  consideration  that  the  metal  strips  are  finite. 

The  use  of  the  effective  dielectric  constant  method  is  not  suffi¬ 
cient  since  it  only  considers  the  propagation  constant  of  the  struc¬ 
ture  without  metal  strips.  From  the  experiments  performed  to  inves¬ 
tigate  the  effect  of  varying  the  y-dlmenslons,  it  seems  difficult 
to  correlate  the  results  for  the  same  structures  having  different 
y-dimenslons.  This  suggests  that  there  exists  large  difference 
between  structures  finite  and  infinite  in  the  y-direction. 


6.2 


The  Experimental  Investigation 


The  apparatus  used  for  the  first  investigation  was  quite 
adequate.  The  waveguide  ref lectometer  system  was  very  effective  and 
the  overall  characteristics  of  the  system  are  quite  stable..  A  piece 
of  apparatus  that  could  have  really  helped  the  measurements  is  a 
sweeper  that  could  go  from  7  to  16  GHz  in  one  sweep,  instead  of 
having  to  use  several  narrower-band  sweeps  to  cover  the  entire  spec¬ 
trum.  Such  a  sweeper  would  have  reduced  the  experimental  time  by 
almost  one-half. 

The  experimental  measurements  made  in  such  a  system  were 
a  function  of  three  variables,  namely,  frequency,  attenuation  or 
reflection  and  magnetization.  A  three  dimensional  plot  would  be 
quite  convenient  for  these  measurements.  Such  plots  would  make  the 
experimental  results  much  more  presentable  and  easily  interpreted. 
Three  dimensional  plots  could  be  produced  if  computer  graphics  had 
been  available.  In  such  a  case  the  experimental  data  of,  say,  one 
magnetization  would  be  directly  stored  in  memory  on  a  fast  sweep 
rate  and  then  sweep  the  magnetization  at  a  slower  rate  than  the  fre¬ 
quency  sweep,  and  putting  information  in  the  computer  of  the  value  of 
the  magnetic  field  by  using  the  external  output  of  a  Gauss  meter. 

If  such  a  system  were  used  it  would  take  only  a  few  minutes 
to  obtain  detailed  experimental  results,  since  the  magnetic  field 
could  be  varied  in  smaller  discrete  steps.  As  It  was,  with  the 
apparatus  used  to  take  the  data  for  the  structures  presented  in  this 
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article,  about  three  hours  were  required  to  plot  the  characteristics 
of  one  structure, 

6.3  Possible  Improvements  and  Uses  of  Proposed  Structure 

From  the  experimental  results  obtained,  it  is  evident  that 
non-reciprocity  exists.  The  only  problem  that  detracts  from  the 
utility  of  the  device  is  that  up  to  now  the  forward  loss  is  about 
10  dB,  This  forward  loss  can  be  reduced  by  a  number  of  techniques 
such  as  adjusting  the  period  and  strip  width  and  the  dimensions  of 
the  guides. 

In  the  future,  an  investigation  of  the  effects  of  period 
and  strip  width  on  the  non-reciprocity  could  be  made  in  much  more 
detail.  One  might  also  consider  using  tapered  metal  strips  to 
reduce  the  reflections  due  to  the  strips  in  the  edge  of  the  device. 

As  long  as  the  non-reciprocity  occurs  in  the  negative 
permeability  region,  the  only  control  on  the  bandwidth  of  the  iso¬ 
lator  is  obtained  by  using  ferrite  materials  with  different  satura¬ 
tion  magnetizations.  This  is  somewhat  limited  since  the  highest 
saturation  magnetization  available  is  in  the  5  K  Gauss  region,  which 
would  produce  a  bandwidth  of  about  3  to  4  GHz. 

In  the  design  of  an  operational  device  one  must  also  take 
into  consideration  the  sizes  of  the  dielectric  and  the  ferrite 
slabs.  In  this  Investigation  most  of  the  work  has  been  made  with 
the  same  thickness  of  ferrite  and  dielectric,  and  the  dielectric 
used  had  a  relatively  high  dielectric  constant. 
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APPENDIX 


1  ...  PROGKAfI  FERLOAIi  (  INPUT  ,  OUT  PUT  ,  PLOTR  ) 

2  ...  C  FERRITE  UAUEGUIDE  UITH  PERIOHIC  METAL  LOADING  IS  ANALTZEU 

3  ...  C  BY  USING  THE  SPECTRAL  DOMAIN  APPROACH. 

A  ...  C  SUBROUTINEISPACE)  CALCULATES  ALL  SPACE  HARMONICS 

5  ...  C  OF  THE  BASIS  FUNCTIONS. 

6  ...  C  SUBROUTINE(UNIFQRM)  SOLUES  UNIFORM  FERR I TE - U I "LECTR I C 

?  ...  C  UAUEGUIDE  UITHOUT  METAL  STRIPS,  AND  ITS  ROOT  IS  USED 

B  ...  C  AS  AN  INITIAL  GUESS  FOR  THE  PERIODIC  STRUCTURE. 

9  ...  C  SUBROUTINETPERIOD)  SOLVES  PERIODIC  FERR I TE- D lELECTR IC 

10  ...  C  WAVEGUIDE  UITH  METAL  STRIPS,  BY  MEANS  OF  MULLER  METHOD. 

11  ...  C  SUBROUTINE(DETER)  CALCULATES  THE  MATRIX  ELEMENTS  OF 

12  ...  C  CHARACTERISTIC  EQUATION  BY  SUMMATION  AND  THEN 

13  ...  C  CALCULATES  ITS  DETERMINANT. 

14  ...  C  SUBROUTINE(IMPED)  FORMALLY  CALCULATES  INPUT  IMPEDANCES/ 

15  ...  C  ADMITTANCES  OF  THE  +-  N-TH  ORDER  SPACE  HARMONICS. 

16  ...  C  SUBROUTINEIEQUIV)  ACTUALLY  CALCULATES  EQUIVALENT 

17  ...  C  transverse  CIRCUIT,  THAT  IS.  F-MATRICES 

18  ...  C  FOR  BOTH  FERRITE  AND  DIELECTRIC  SLABS  AND  THE 

19  ...  C  CHARACTERISTIC  ADMITTANCE  FOR  AIR  REGION  FOR  ANY  BETA, 

20  ...  C  AND  THEN  OBTAINS  INPUT  IMPEDANCE  FOR  METAL  SURFACE. 

21  ...  C  SUBROUTINE  (GRAPH)  PLOTS  THE  VALUES  OF  KD/PAI 

22  ...  C  VS.  BD/PAI  AND  AD/PAI 

23  ...  C 

24  ...  DIMENSION  C 1  (  50 ) , C2 ( 50 ) , C3 ( 50 ) , C4 ( 50  ) 

25  ...  DIMENSION  AA(50)  ,BB<2, 50) ,CC(2, 50) ,DD(2, 50) .EE(6, 50) 

26  ...  COMPLEX  SETA, ERROR, BP, BK,HILU,HI2, AM, AK, EM, SI 

27  ...  PAI=3. 1415926 

28  ...  C  WAVEGUIDE  STRUCTURE  IS  GIVEN  HERE. 

29  ...  C  NOTE  STRUCTURE  13  INFINITE  IN  THE  Y-DIRECTION 

30  ...  C  EF  =  DIELECTRIC  CONSTANT  OF  FERRITE  LAYER 

31  ...  C  HI  =  DC  MAGNETIC  FIELD  APPLIED  IN  Y-DIRECTION  (K  OERSTEDS) 

32  ...  C  SM  »  SATURATION  MAGNETISATION  OF  FERRITE  ( K  GAUSS) 

33  ...  C  TF  =  THICKNESS  OF  FERRITE  LAYER  (MM) 

34  ...  C  ED  =  DIELECTRIC  CONSTANT  OF  DIELECTRIC  LAYER 

35  ...  C  TD  =  THICKNESS  OF  DIELECTRIC  LAYER  (MM) 

36  ...  C  D  =  PERIOD  OF  METAL  STRIPS  (MM) 

37  ...  C  A  =  WIDTH  OF  METAL  STRIPS  (MM) 

38  ...  C  UL  =  LINE  WIDTH  OF  FERRITE  LAYER  (K  OERSTEDS) 

39  ...  READ  100,EF,HI,SM,TF,ED,TD,D,A.UL 

40  . . .  100  F0RMAT(9F8) 

41  ...  C  PRINTING  OF  ALL  THE  DATA 

42  ...  PRINT  150 

43  ...  1  50  FORMAT!  1  HI  , 'FERRITE  - '/4X  ,  PERM  I  T I V  I  T  Y  ,  4  X  , 

44  ...  1'INT.  MAG.  FIELD  (KOE)' 

45  ...  1,4X,  SATURATION  (KGAUSS)  , 4X ,  IH I CKNESS  (MM)', 

46  ...  15X,  LINE  WIDTH  (KOE)') 

47  ...  PRINT  200,EF,HI,SM,TF,UL 

48  ...  200  FORMATdH  ,  4X  ,  F7 . 3 , 1  3X  ,  F7 . 3  ,  1  7X  ,  F7 . 3  ,  1  3X  ,  F7 . 3 , 1  3X  ,  F  7 . 3  ) 

49  ...  PRINT  250 

50  ...  250  FORMAT ( t HO,  DIELECTRIC  - ' / 4 X , ' PERM  I T I V I T Y  '  , 4 X , 

51  . . .  1  THICKNESS  ( MM ) ' ) 

52  ...  PRINT  275,ED,TD 

53  ...  275  FORMATdH  ,  5X  ,  F  7 . 3 , 8X  ,  F7 . 3  ) 

54  ...  PRINT  800 

55  ...  800  FORMATdHO,  METAL  STRIPS  - ’ / 4X , ' PER lOD  (MM)', 

56  ...  I 4X ,  UI DTH  (MM )  ) 

57  _  PRINT  850, D, A 

58  ...  850  FORMAT((H  , 4X , F7 . 3 , BX , F7 . 3 ) 

59  ...  DA«D/A 

60  ...  C  NO.  OF  SPACE  HARMONICS 
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60  ...  C  NO.  OF  SPACE  HARMONICS 

61  ...  C  NMAX  MUST  BE  LESS  THAN  OR  EQUAL  TO  50. 

62  ...  NMAX=40 

63  ...  C  FOURIER  TRANSFORM  OF  BASIS  CURRENT  DISTRIBUTIONS 

64  ...  C  UILL  BE  CALCULATED. 

65  ...  C  BASIS  FUNCTIONS  ARE  INDEPENDENT  OF  FREQUENCY  AND  UAOEGUIDE 

66  ...  C  STRUCTURES  AND  DETERMINED  BY  ONLY  PERLPD  AND  STRIP  UIDTH. 

67  ...  CALL  SPACE(DA,C1 ,C2,C3,C4,NMAX) 

68  ...  C  NO.  OF  BASIS  FUNCTIONS 

69  ...  C  N  MUST  BE  LESS  THAN  OR  EQUAL  TO  4. 

70  ...  N=4 

71  ...  C  CONOERGENCE  CRITERIA 

72  ...  EPS=1.E-9 

73  ...  PRINT  400, N, EPS 

74  ...  400  FORMATdH  ,  N0.  OF  BASIS  FUNCTIONS  ='.I2/X, 

75  ...  ♦'CONVERGENCE  CRITERIA  =',E13.5//) 

76  ...  C  2.8  IS  GYROMAGNETIC  RATIO.  (UNIVERSAL  CONSTANT) 

77  ...  UL=UL/2. 

78  ...  HILU  =  CMPLX(HI ,UL) *  2.8 

79  ...  SM=SM^2.8 

80  ...  HI2=HILUfHILU 

81  ...  C  LOOPING  ON  THE  +VE  AND  -VE  DIRECTIONS  OF  THE  MAGNETIC  FIELD 

82  ...  DO  10  1=1,2 

83  ...  C  THE  DIRECTION  OF  MAGNETIZATION  IS  SPECIFIED. 

84  ...  HI=-HI 

85  ...  HILU=-HILU 

86  ...  SM=-SM 

87  . . .  IF(1 .EQ. 1 )  GO  TO  1 

88  ...  PRINT  110 

89  ...  110  FORHATdHO,  POSITIVE  DIRECTION  /) 

90  . . .  GO  TO  2 

91  ...  1  PRINT  120 

92  ...  1  20  FORMATdHO,  NEGATIVE  DIRECTION  /) 

93  ...  2  CONTINUE 

94  ...  C  FREQUENCY  UILL  BE  SWEPT. 

95  .  .  .  DO  20  J  =  6, 18 

96  ...  K=J-5 

97  ...  F=FLOAT(J) 

98  ...  F2=F»F 

99  ...  C  TENSOR  PERMEABILITY  OF  MAGNETIZED  FERRITE  IS  GOVERNED 

100  ...  C  BY  POLDER  S  FORMULA,  AS  SHOWN  BELOW. 

101  ..!  C  AM  IS  DIAGONAL  ELEMENT  (MYU)  . 

102  ...  C  AK  IS  OFF-DIAG.  ELEMENT  (KAPA)  . 

103  ...  AM=1 .♦HILU»SM/<HI2'F2) 

104  ...  AK=SN^F/<HI2- F2) 

105  ...  C  PERMIABILITY  TENSOR  BELOW  SATURATION 

106  ...  C  EQUATIONS  DUE  TO  GREEN  AND  SANDY 

107  ...  SMX  SMd75/2.8 

108  ...  IF < ABSTHI ) .GE.ABS(SN) )  GOTO  111 

109  ...  IF < ABS(HI ) .GT . ABS(SMX) )  GOTO  111 

110...  PRINT  9<J9  ,F 

111  ...  999  FORMAT  (IH  ,F15.7) 

112  ...  UM0  =  2./3.^SQR1  1 .-(SM/F) •  ♦Zitl  ./3. 

113  ...  AM  =  UM0*(1.-UN0  )  ‘SORT  dHI»2.87SM)K3) 

114  ...  AK=-HI*2.8/F 

115  ...  Ill  CONTINUE 

116  ...  C  EM  IS  EFFECTIVE  PERMEABILITY 

117  ...  C  SI=)!AFA/MYU 

118  ...  C  UN  IS  WAVENUMBER  IN  FREE  SPACE. 

119  ...  EM=<AM»AM-AK»AK)/AM 

120  ...  SI=AK/AM 

121  ...  UN=F^2.^PAI/300. 

122  ...  PRINT  300  ,F,EM,SI 

1  23  ...  300  FORMATdH  ,  '  F  =  •  ,  F  7 . 3 ,  '  t  GH  Z  >  '  ,  2X  ,  EM=  ,  E  1  3 . 5 , 2X  ,  E  1  3 . 5 , 2X  , 

124  ...  t'SI=',EI3.5.2X,E13.5) 

125  ...  PRINT  350,HILU.SM,AM,Ai; 
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1  25 
1  26 
1  22 
1  28 
1  29 
1  30 
131 
1  32 
1  33 
1  36 
1  35 
1  36 
1  37 
138 
1  39 
1  40 
141 
1  42 
I  43 
1  44 
P45 
1  46 
I  42 
1  48 
149 
1  50 

151 

152 

153 

154 

155 

156 
1  57 
1  58 

159 

160 
161 
162 
1  63 
164 
1  65 
1  66 
1  67 
1  68 
1  69 
1  70 
171 
1  72 
1  73 
1  7  4 
1  75 
1  76 
1  7  7 
1  78 
1  79 
1  80 
181 
1  82 

183 

184 
1  85 
1  86 
1  87 
1  38 
1  8° 
1  90 


PRINT  350,HILU,SM,AM,AK 

350  FORMATdH  , -"H 1  LU=  ,  P?  •  3 , 2X  ,  F7 . 3 , 2X  ,  '  SM^  '  ,  F7 . 3  ,  / 

1  ,  1  X,  AM  =  '  ,F7.3,2X,F7.3,2X,  'Ai;=  ,F7.3,2X,F7.3) 

Uli  =  UNi  D/PAI 
PRINT  351. UD 

351  FORMATdH  , 'KD/P  A I  =  '  ,  E  1  3 . 5  ) 

C  AT  FIRST,  UNIFORM  UAVEGUIDE  IS  ANALYZED. 

CALL  UNIFQRM(EF,EH.3I,TF,ED,Tli,WN,ri,IiO) 

C  BO  IS  USED  FOR  INITIAL  GUESS  FOR  PERIODIC  STRUCTURE. 
BOK=BO/UN 
PRINT  355, BOK 

355  FORMATdH  ,  '  BO/UN=  '  ,  E  1  3 . 5  ) 

C  PERIODIC  UAVEGUIDE  UILL  BE  SOLVED. 

CALL  PERI0D<UN,EF,EM,SI,TF,ED,TD,D,N,C1 ,C2,C3,C4, 
*NMAX,BO,BETA,MI ,EPS,ERROR ) 

C  BETA  IS  PROPAGATION  CONSTANT  OF  PERIODIC  UAVE . 

C  REAL  PART  IS  PHASE  CONSTANT 
C  IMAG.  PART  IS  ATTENUATION  CONSTANT  DUE 
C  TO  REFLECTION  OR  LEAKAGE. 

C  MI  IS  NO.  OF  ITERATION  IN  MULLER  METHOD. 

C  ERROR  IS  DEFINED  BY  FUNC T I  ON ( BE T A ) . 

PRINT  360. BETA. MI, ERROR 

360  FORMATdH  ,  '  BET  A=  ' ,  E  1  3 . 5 , 2X  .  E  1  3 . 5 , 5X  ,  M I  =  '  ,  1  3  , 

♦  2X ,  ERRORS  ,E1 3.5. 2X, El  3.5) 

BK=PETA/UN 
PRINT  370, BK 

370  FORMATdH  ,  '  BET  A/UN=  ,  E  1  3 . 5 , 2X  ,  E  1  3 . 5  > 

BP  =  BK  »UD 

380  FORMATdH  ,  '  BD/P  A I  =  •" ,  E 1  3 . 5 . 2X  ,  E  1  3 . 5/ ) 

PRINT  380, BP 

C  FORMATION  OF  ARRAYS  TO  MAKE  PLOYING  EASYER 
BD1 =BO*D/PAI 
BD2=REAL(BP) 

BD3=ABS<AIMAG(BP) ) 

C  LOGIC  TO  AVOID  GOING  OUT  OF  RANGE  ON  THE  PAPER  UHEN  PLOTING 
IF  (UD.LT.O.O)  UD=0.0 
IF  (UD.GT.1.5)  0D=1.5 
IF  (BD1.LI.0.0)  BD1=0.0 
IF  <BD1.GT.2.5)  BD1=2.5 
IF  ( BD2.LT.0.0 )  BD2=0. 0 
IF  <BD2.GT.2.5)  BD2=2.5 
IF  (BD3.lt. 0.0)  PD3-0.0 
IF  (BD3.GT.2.5)  Bti3  =  2.5 
AA(K)-UD 
BB( I ,K>=BD1 
CC(I,K)=BD2 
DD( I ,K)=BD3 
20  CONTINUE 
10  CONTINUE 
DO  15  1=1,2 
DO  25  K=1 ,13 
EEC  I .K)=BB( I ,K) 

EEC  d2,K)=CC(  I  ,K  » 

EE ( 1*4 ,K )=DDC I ,K) 

25  CONTINUE 
15  CONTINUE 
UL-=UL  *2. 

SM=SM/2.8 

CALL  GRAPH  ( EE , AA , EF , H I , SM . TF , ED , TD , D , A , Ut) 

STOP 

END 

C 

C 

c 

c 

c 
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1  90 
191 
1  92 
1  93 
1  9A 
1  95 
1  96 
1  97 
1  98 

1  99 
200 
201 
20  2 

203 

204 

205 

206 

207 

208 

209 

210 
21  1 
212 

213 

214 

215 

216 
217 

2  I  8 

219 

220 
221 
6  *1  *> 

223 

224 

225 

226 
no  7 
228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 
*152 

25  3 

254 

255 


C 

SU8h:0iJ(INE  SPALElfD.CI  ,  C2  ,  C3  ,  C4  .  NMAX  ) 

DIMENSION  Cl (50) ,02(50) ,C3( 50) ,04(50) 

C  SPACE  HARMONICS  OF  FIELD  D I  SI R I  DOT  I ONS  ON  THE  APERTURE. 
PAI^-3.  1415926 
X  =  PAIK  1.-1  ./PD) 

C  EDGE  CONDITION  IS  NOT  CONSIDERED. 

C  THE  BASIS  FUNCTIONS  USED  ARE 

C  El  =  PAI »D/2/(D-A)*EXP( J»BOtX>iCOS(PAI kX/'D-A) ) 

C  E2  =  J ►PAIHi/(D-A) tEXP( J‘BO> X ) *SIN(2»PAI ♦X/( D-A) ) 

C  E3  ^  -3iPAItD/2/(D-A)*EXP( J»b0iX)*C0S(3tPAI IX/(D-A) ) 

C  E4  -  L’kpAI  tDtJ/(D-A)«EXP(J*B0»X)»SIN(4tPAIi  X/'(D-A)  ) 

C 

C  THE  SPACE  HARMONICS  OF  THE  BASIS  FUNCTIONS  HAVE  BEEN 
C  CALCULATED  AN AL Y I T  I C ALL Y  AND  ARE  GIVEN  BY  C 1  < N ) , C 2 ( N ) , 

C  C3(N).C4(N)  RE5FECTIVLY 
DO  10  N=1 .NMAX 
XN  =  X:tPLOAT(N) 

XN2=XN*XN 
E  =  XN2M./PAI/PA1 
A=1 .-E 
B=1 .-E/4. 

C=1  .  -E/9. 

D=1 .-E/t6. 

SN=SIN(XN) 

CN=CnS(XN) 

Cl (N)=CN/A 
C2(N)=SN/B 
C3 ( N ) =CN/C 
C4(N)=SN/D 
10  CONTINUE 

PRINT  100, NMAX 

100  FORMAKIH  ,  NO .  OF  SPACE  HARMONICS  /X,  NMAX^  ,12) 

RETURN 

END 

C 

C 

c 

c 

c 

SUBROUTINI  PERI0Ii(UN,EF.EM.SI,TF,ED.TD,E,N,C1.C2,C3.C4, 
»NM AX. BO, BET A, MI, EFT, ERROR) 

COMPLEX  BETA. B(3),F(3), ERROR. F32,F2I ,) 321 .LAMDA, 
»DET.RDET,TP,  TN,II.BN,EM,SI 
DIMENSION  C1(50).C2(50),C3(50),C4(50' 

C  PERIODIC  STRUCTURE  IS  SOLVED. 

EPS^-EPT 
P  =  E 
MC  =  1 

C  INITIAL  GUESS  OF  BETA 
B( 1 )=B0 

B( 2 ) ^B( 1 ) »0.95 
B(3)=B( 1 )t0.9 

C  PROPAGATION  CONSTANI  OF  EIGEN  NODE  UILL  BE  SOLVED 
C  BY  MEANS  OF  MULLER  METHOD. 

C  MULLER  S  METHOD  IS  A  GOOD  METHOD  TO  FIND  THE  ROOTS  OF 

C  A  COMPLEX  FUNCTION  IF  THE  DERIVATIVE  OF  THE  FUNCTION 

C  IS  NOT  KNUUN.  FOR  MORE  DETAIL  ABOUT  MULLER'S  METHOD 
C  SEE  -ELEMENTS  OF  NUMERICAL  ANALYSIS'  BY  PETER  HENRICI, 

C  JOHN  UILEY  a  SONS.  1964.  CHAPTER  10  ,PAGE  198. 

C  OBJECTIVE  FUNCTION  IS  APPROXIMATED  BY  POLYNOMINAL 

C  OF  SECOND  DEGREE. 

DO  10  1=1,3 

C  (DEIER)  IS  A  SUBROUTINE  WHERE  YOU  FEED  IN  BETA  AND  IT 
C  UILL  GIVE  YOU  BACK  THE  DETERMINANT  F,UHICH  WE  NEED  TO 
C  MINIMIZE  TO  ZERO  IF  POSSIBLE. 
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2tiii 

L 

niNI«12E  TO  ZERO  IF  POSSIBLE. 

2S4 

C 

CALLING  OF  DETER  FOR  THE  FIRST  THREE  INITIAL  POINTS 

257 

C 

OF  HOLLER  S  HETHQD 

258 

CALL  DETER  <B(  I  I.FUI.EF,  EH,  SI,  TF,  ED,  TD,  ON,  P.N, 

25V 

1C1 ,C2.C3,C4,NHAX) 

260 

10  CONTINUE 

261 

30  IF (CABS (B<3)-B(2) ) .LE.EPS)  GO  TO  24 

262 

F32=(F<3)-F(2)>/(B(3)-B(2)> 

263 

IF < CABS ( B ( 2)-B( 1 >) .LE.EPS)  GO  TO  24 

264 

F21=(F(2)-F(I))/(B(2)-B(1)) 

265 

IF  (  CABS  (  B<  3)-B<  t  ))  .LE.EPS)  GO  TO  24 

266 

F321=<F32-F21)/(B(3)-B(1)) 

267 

LAHEiA=(B(3)-B(2))»F321*F32 

268 

IF(CABS(LAHDA ) .GE. 1 .E37 )  GO  TO  27 

269 

DET=LAHDA»*2-4.»F(3)*F321 

270 

RDET=CSQRT(DET) 

271 

TP=LAHDA+RDET 

»)  ■»  n 

TN=LAHDA-RDET 

273 

IF(CABS(TP) .GE.CABS(TN) )  GO  TO  22 

274 

H=-2.»F<3)/TN 

275 

GO  TO  23 

276 

22  CONTINUE 

277 

H=-2.iF<3)/TP 

278 

23  BN=B(3)tH 

279 

C 

NEU  APPKOXIHATE  ROOT 

280 

B  (1  ) =  B  <  2 ) 

281 

B(2'=B(3) 

282 

B(3)=BN 

283 

F< 1 )=F(2) 

284 

F(2)=F<3) 

285 

IF<CABS<H/B<3)).LE.EPSJ  GO  TO  40 

286 

IF(MC.G£.50)  GO  TO  24 

287 

HC=HC»1 

288 

C 

LOOPING  ON  DETER  TO  OPTIHIZE  THE  SOLUTION  OF  HAOEING  F=0 

289 

CALL  DETER<B(3> ,F<3) ,EF,EH,SI ,TF,ED, rD,UN,P.N. 

290 

1C1 ,C2.C3.C4,NrtAX) 

291 

GO  TO  30 

292 

24  HI=HC+50 

293 

GO  TO  50 

294 

40  HI=nC 

295 

50  BETA=B(3) 

296 

C 

FINAL  CALL  OF  DETE  TO  PRINT  RESULT 

n  tf-j 

CALL  DETER<B(3),F<3),EF,EH,SI,TF,ED,TD,UN,P,N, 

293 

1C1 ,C2.C3,C4,NMAX) 

2?  y 

ERR0R=F(3) 

300 

RETURN 

301 

27  MI=R99 

302 

ERRUR-I .E49 

302 

RETURN 

304 

END 

305 

C 

306 

C 

307 

C 

30H 

C 

309 

C 

310 

SUBROUTINE  DETER < B , F , EF , EH . S I , 1 F . ED , TD , UN , R , H , 

3  1  1 

1C1 ,C2,C3.C4,NMAX) 

312 

C 

CALCULATION  OF  DETERHINANT  OF  CHARACTERISTIC  EQUATION 

313 

D I  HENS  I  ON  Cl (50) ,C2<50) .C3(50 ) .C4(50 ) 

314 

COHPLEX  B,F.PI1,P12.P13.PI4,P22,P23,P33,P24,P34,P44 

315 

»,D(4).Z,2P,ZN,ZA,ZS,0.EH,SI 

316 

C 

CALCULATION  OF  HATRIX  ELEHENTS 

317 

a=R 

318 

C 

(IHPED)  GIVES  THE  VALUES  OF  THE  ADMITANCES  AT  THE  POINTS 

319 

c 

UHERE  THE  HETAL  STRIPS  ARE  SITUATED  THAT  IS  YU*YL  . 

320 

c 

NOU  ZP  AND  PN  ARE  THE  ADHITANCES  OF  THE  fVE  AND  -VF 
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320 

C 

NUU  Zf  AND  PN  ARE  THE  ADMITANCES  OF  THE  +0E  AND  -VE 

321 

C 

DIRECTIONS  RESF-ECTIVLY. 

322 

c 

SO  FINALLY  UE  UILL  END  UP  UITH  2fNMAXH  SPACE  HRMONICS. 

32  3 

c 

NOTE  THAT  THE  MATRIX  IS  FORMED  PY  THE  PRODUCT  OF 

32‘1 

c 

EN'ENiYN  AND  NOT  BY  JN*JN*ZN. 

325 

c 

326 

c 

THE  ZEROTH  SPACE  HARMINIC  IS  CALLED  HERE 

327 

CALL  IMPED(ZP,ZN.B,0,EF,EM.SI . TF , E D , TD , UN , 0 ) 

328 

Z  =  ZP 

329 

P 1  1  =  Z 

330 

PI  2= <0.0, 0.0) 

331 

PI  3=^2 

332 

P11=<0. 0.0.0) 

333 

P22=(0. 0,0.0) 

331 

P23=(0.0.0.0) 

335 

P21-  <0.0.0.0) 

336 

P33  =  Z 

337 

P31-(0. 0,0.0) 

338 

P41=(0. 0.0.0) 

339 

DO  10  N-  1  ,Ni1AX 

310 

c 

LOOPING  ON  THE  REST  OF  THE  SFACL  HARMONICS 

311 

c 

TO  FORM  THE  MATRIX. 

312 

CALL  IMPED(ZP,ZN,B.N,EE,EM,SI  ,  TF  ,  ED  .  T  D  ,  U  N  ,  t) ) 

313 

ZA=ZP*ZN 

311 

ZS=ZP-ZN 

315 

D  1  -  C  1  (  N  ) 

316 

D2=C2<N) 

317 

D3^C3<N) 

318 

D  1  =  C  4  ■  N  ' 

319 

PI  1  -^Pl  1  pZApDI  »Dt 

350 

PI  2-PI  2  +  ZSHil  •  D2 

351 

PI  3  =  PI 3  »iA*D1 ■ D3 

352 

PI  1=P1 1  +  ZS*D1 fDI 

353 

P22=  P22tZA»D2  Hi2 

351 

P23-P23»ZS>D2tD3 

355 

P21^P21+ZA*D2tD1 

356 

P33-P33+ZAID3ID3 

357 

P31  =  F  31t2S*D3tD1 

358 

P11  =  P1H-ZA*D1»D1 

359 

10  CONTINUE 

360 

c 

CALCULATION  OF  DETERilINANT 

361 

D< 1 )-P1 1 

362 

D(  2) =P1 1 7P22-P1 2'r  i 2 

363 

D(  3  )  =  F'1  1  7P22«F  33* 2.  fP1  2*F23<P1  3-PI  3*P1  3*P22-P23< 

364 

tP23*F-1  1  -PI  2fP1  2*P33 

365 

IK  1 )  --P1  1  MP22>  (F'33*P11-P31  fP31 )  -  P23t  (  P23  t  P11-P31  • 

366 

1  P21  )  ‘P211  (  P23  •  P51-P33  fP21 )  )  -  F  12  1  (  PI  2  M  P33  •  P11-F' JH 

367 

2P31  )-P23  *  (PI  3<P11-P31  'PI  1 )  +  P21  »  (PI  3»P  31-P331P1  -1  )  )♦ 

368 

3P1  3  KP1  2  ■  <P23*  FH-P  51  'F  -'l  )-P22K  P  i  ;■  tP41  -  P31  tP1  1  > 

369 

1Pf21KF13tP21-P2  3*P11))-P11KP12KP23Pp31-P3;HI'21)- 

370 

5P22' (PI  3  >P31-P33»P1 1 )+P23t (PI 3»  P21-P231F1  1 ) ) 

371 

F  =  D  <  M  ) 

672 

return 

373 

END 

371 

c 

375 

c 

376 

c 

3  77 

c 

3  7  8- 

r 

3  79 

SUPR8UI  INE  IMPFIKZP.ZN.B.N.LF.LN.SI  .  Tl  .  ED  .  T  D  .  UN  ,  P  ) 

380 

c 

EOALUATION  Of  INfUI  IMPEDANCE  Of  RfACE  HnKMUNICS  BY  USING 

331 

c 

EOUIOALENT  TRANSMISSION  LINL. 

382 

COMPI  EX  B.BNP.BNN,ZP.ZN,Z.EM,SI 

183 

PA  I ^3. 1 115926 

381 

BNP  B»2. >f  AI, f  *FLOAT(N) 

385 

BNH=^  H  -  2  .  1 1  A  I  'p  •  FLOAI  '  N  ) 

120 


385 

PNN  =  li-2.»PAX/l-<FL0AT<N) 

386 

C 

INPUT  IHHITANCE  FOR  +N-TH  ORDER  HARMONICii. 

387 

CALL  EQUIV<BNP,EF,EM,SI,TF,ED,T[i,UN,2  ) 

388 

ZP  =  Z 

389 

C 

INPUT  IMHITANCE  FOR  -N-TH  ORDER  HARMONICS. 

390 

CALL  EQUIV(BNN,EF,EM,SI ,TF ,ED,TD,UN,Z) 

391 

ZN  =  2 

392 

RETURN 

393 

END 

39-) 

C 

» 

395 

C 

396 

C 

397 

C 

398 

C 

399 

SUBROUTINE  ECU  I  V ( B , EF , EM , S 1 , TF , ED , T D . UM  ,Z) 

400 

C 

INPUT  IMMITANCE  OF  SPACE  HARMONICS  IS  CALCULATED. 

401 

C 

AIR-FERRITE  -  METAL  DIELECTRIC  -AIR  STRUCTURE 

402 

COMPLEX  B,B2,PF.PA.PP.SF,SG,AF,BF,CF,DF,AD,BD.C 

403 

1YL.Y,YF,YD,Z,U,SD.YU.AU.BU,CU,DU,EM,SI 

404 

U=(0. . 1 . ) 

405 

B2=B*B 

406 

UN  =  UM 

407 

UN2--UN*UN 

408 

C 

SELECTION  OF  BRANCH 

400 

PA^CSC)RT(B2-UN2)/U 

4  1  0 

PF=CSQRT<UN2'EFiEM-B2) 

4  1  1 

FP-CSQRT(UN2*ED-B2> 

412 

SF=PF»TF 

413 

SD=PD*TP 

4  1  4 

SG-^SI^B/PF 

415 

YF=PF/gN/EM 

4  1  6 

YD=PD/UN 

417 

C 

F-MATRIX  OF  FERRITE  SLAB 

4  1  8 

AF-CC0S(SF)-S6*CSIN<SF ) 

410 

DF=CC0S(SF)+SG*CS1N(SF> 

420 

BF^U'CSIN(SF)/YF 

421 

CF=U*CSIN(SF)»YF*(1 .+SGkSG) 

4  22 

c 

F-MAfRIX  OF  DIELECTRIC  SLAB 

423 

AD=CCOS(SD) 

424 

DD  =  AIi 

425 

BD=U*CSIN(SD)/YD 

426 

CD-U*CSIN(SD) < YD 

4  2  7 

c 

CHARACTERISTIC  ADMITTANCE  OF  AIR  REGION 

428 

Y=PA/UN 

429 

c 

INPUT  ADMITTANCE  OF  LOUER  REGION 

430 

YL=(CB+DDtY)/<AD+BD*Y) 

431 

c 

INPUT  ADMITTANCE  OF  UPPER  REGION 

432 

YU=(CF*AF»Y)/(DF*BF*Y) 

433 

c 

INPUT  IMPEDANCE 

434 

Z=YU+YL 

435 

RE  TURN 

436 

END 

437 

c 

438 

c 

439 

c 

440 

c 

441 

c 

442 

SUBROUTINE  UNIFORri  (  EF  ,  EM .  S  F,  TF ,  ED  ,  TD  ,  UN .  D, BO  > 

443 

c 

DISPERSION  RELATION  OF  UNIFORM  FERR I TE - D I ELECTR IC 

444 

c 

UAVEGUIDE  UITHOUT  METAL. 

445 

COMPLEX  B,Z.Y,£M,SI 

446 

PAI»3. 1415926 

447 

c 

determination  of  the  region  for  THE  SEARCH 

448 

c 

OF  A  SOLUriON  FOR  BETA. 

449 

BM=EFtREAL<EM) 

450 

BN=EDt3. 

121 


450 

451 

452 

453 

454 

455 
454 
45? 

458 

459 

440 

441 

442 
4  4  3 

444 

445 
444 
44? 

448 

449 
4?0 
4?1 
4  .'2 
4?3 
4?4 
4?5 
4?4 
4?? 
4?8 
4?9 

480 

481 

482 

483 

484 

485 
484 
48? 
4  18 

489 

490 

491 

492 

493 

494 

495 
494 

4  9? 

498 

499 

500 

501 

502 

503 

504 

505 
504 
50? 
508 

50  9 
510 

51  1 
512 

5  I  3 

514 

515 


BN=ED»3. 

IFv'BM.Lr.BN)  I(M  =  ttN 
EiM  =  S0RT<BM)  ►0.99999»J. 

DB=(BM-1 .00001 1/100. 

ZI=0.0 
BD  =  0.0 

DO  30  K^1 ,101 

B=(  1  .  +  BBtFL0i'1T(K-l  1  I  *UN/3.0 

CALL  EOUIOIB.EF.LM.SI.TF.EIr,  IIi,UN,n 

Z  =  Y 

IFIK.EO.  I  1  GO  TO  1 
4  ZJ--AIMAG(Z) 

BP  =  REAL  <  B 1 
S=Z I fZ J 

IFIG.LT.O.O)  GO  TO  3 

ZI=ZJ 

BR-BP 

GO  TO  2 

1  ZI=AIhAG<Z) 

BK=REAL<B> 

GO  TO  4 

3  IF ( OBSIZI ) .GT. 1 . )  GO  TO  4 

C  INTERPOLATING  FOR  A  BETTER  VALUE  OF  BETA 
C  ONCE  THE  SOLUTION  CHANGES  FROM  +VE  TO-VE 
BIi=(ZI»BP-ZJtBR)  Hi/(ZI-2J)/PAI 
PRINT  300, SD 

300  FORMATCIH  .  BD/PAI=  ,E13.5) 

4  ZI=ZJ 

BR  =  BP 

2  CONTINUE 
30  CONTINUE 

BO  =  tIi*PAI/D 
PRINT  200, BO 

200  FORMATdH  , BO*-' ,  E  1  3 . 5/X  ,  UN  I F  ORM  END'.' 

RETURN 

END 

C 

C 

C 

C 

C 

SUBROUTINE  GRAPH  ( EE , AA , EF , H I , SM , TF , ED , TD , D , A . UL 1 
C  (GRAPH)  PLOTS  THE  VALUES  OF  KD/PAI  VS.  BD/PAl  AND  AD/PAI 
C  THAT  IS  THE  REAL  AND  IMAGINARY  PARTS  OF  BE i A 
C  FOR  THE  FERRITE-METAL  STR I  PS- D lELECTR I C  STRUCTURE. 

C  IT  ALSO  PLOTS  THE  SOLU1ION  OF  BETA  FOR  JUST  THE  FERRITE- 
C  DIELECTRIC  STRUCTURE  THAT  IS  UITH  NO  METAL  STRIPS. 
DIMENSION  EE(4.50).AA(50).EE1(50) 

CALL  PLOTS  (0,0,5LFL0TR) 

C  PRINTING  OF  THE  DATA 

CALL  SYMBOL  ( 2 . 0 , . 4 . . 1 , 2 , 0 . . -  1  ) 

CALL  SYMBOL (2.2, .55, . 1 . 

♦  33H-  BETAHi/PI  UITH  NO  METAL  S  TR  IPS  ,  .  0 , 3  3  1 
CALL  SYMBOL  ( 2 . , . 35 , . 1 , 3 . 0 . . - 1 > 

CALL  SYMBOL  (2.2,.3,.1, 

»39H-  BETAHi/PI  UITH  PERIOPIC  METAL  STR I  PS , 0 . , 39 1 
CALL  SYMBOL  ( 2 . , . 1 , . 1 , 5 , 0 . . -  1  1 
CALL  SYMBOL  (2.2,.05,.1, 

t40H-  ALPHA»D/PI  UITH  PERIODIC  METAL  STR I  PS , 0 . , 4 0 1 
CALL  SYMBOL  (  1  .  5 , 5  .  ?5 ,  .  1  ,  1  I  HE  PS  I  L  ON- F  =,0.,in 
CALL  NUMBER  ( 2 . 4 , 5 . 75 , . 1 , EF , 0 . , 2 ) 

CALL  SYMBOL  ( 4 . 5 , 5 . ?5 , . 1 , 1 OHD  =  MM,0.,10) 

CALL  NUMBER  ( 4 . 8 , 5 . 75 , . 1 , U , 0 . , 2 1 

CALL  SYMBOL  ( 4 . 5 , 5 . 75 , . 1 , 1 1 HTF  =  MM.0.,11) 

CALL  NUMBER  ( 4 . 9 , 5 . 75 , . 1 , TF , 0 . . 2 ) 

CALL  SYMBOL  ( 1 . 5 . 5 . 5 , . 1 . 11  HEPS  I  LON -D  =,0.,11) 


122 


515 

516 

517 

518 

519 

520 

521 

522 

523 

524 


CALL  SYH60L  (  1 . 5 , 5 . 5  ,  .  1  ,  1  1  HEPS  I LON-D  ,  0  .  ,  1  1  ) 
CALL  NUHBER  ( 2 . 6 , 5 . 5 , . 1 , E U , 0 . , 2 ) 

CALL  SYMBOL  ( 4 . 5 , 5 . 5 , . 1 , 1 OHA  =  MM,0.,10) 

CALL  NUMBER  ( 4 . 8 , 5 . 5 , . 1 . A , 0 . , 2 ) 

CALL  SYMBOL  < 6 . 5 , 5 . 5 , . 1 , 1 1 HT D  =  MM, 0., Ill 

CALL  NUMBER  ( 6 . 9 , 5 . 5 , . 1 , T D , 0 . . 2 > 

CALL  SYMBOL  ( 1 . 5 , 5 . 25 , . 1 , 

►2eHBC  MAGNEriC  FIELD  =  K  0E,0..28) 

CALL  NUMBER  ( 3 . 4 . 5 . 25 . . 1 , H I , 0 . , 2 ) 

CALL  SYMBOL  (5.5,5.25,.!, 


J  4.  O 

526 

527 

528 

529 

530 

531 

532 

533 

534 

535 

536 

537 

538 

539 

540 

541 

542 

543 

54  4 

545 

546 

547 

548 

549 

550 

55  1 
55  2 


5  6^ 

568 

569 


*26HH-SArURAT ION  =  K  GAUSS, 0. ,26) 

CALL  NUMBER  ( 6 . 9 , 5 . 25 , . 1 , SM . 0 .  ,  2  ) 

CALL  SYMBOL  (1.5.5.0,.l, 
t29HFERRITE  LINE  WIDTH  =  K  0E,0.,29) 

CALL  NUMBER  ( 3 . 5 , 5 . 0 . . 1 , UL , 0 . , 3 ) 

CALL  SYMBOL  ( 6 . 5 , 1 . 75 , . 07 . 

*26H -  NEGATIVE  D I REC T 1 0 N , 0 . , 26 ) 

CALL  SYMBOL  ( 6 . 5 , I . 55 , . 07 . 

.  +268-  -  -  -  POSITIVE  DIRECTION, 0. .26) 

CALL  AXIS  ( 1 .0, I .0, 1 7HBD/P1  AND  AU/PI , -0 , 7 . 5 , 
*0.0,0.0,0.3333333) 

CALL  AXIS  <1 .0,1 .0,5HKD/PI,5.4.5,90.0,0.0,0.3233333) 
CALL  PLOT  ( 1 .0, 1 .0, 3) 

CALL  PLOT  (4.0,4.0,21 
CALL  PLOT  (7. 0,1. 0,2) 

SETTING  OF  THE  ORIGIN 
CALL  ORIGIN  (1.0,1.0,01 
DO  11  N=1,6 
DO  9  K=1 . 13 
EE1 (K)«EE(N,K( 

9  CONTINUE 

SCALING  OF  THE  AXIS  THAT  IS  UNITS  PER  INCH  OF  PLOT 
AA< 1 4 )=0,0 
EE1 <14)=0.0 
AA(  1  5)  =  0. 33333333 
EE1 (151=0,33333333 

LOGIC  TO  DETERMINE  THE  APPROPRIATE  SHADES  AND 
SYMBOLS  OF  GRAPHS. 


553  ... 

IF 

(N.EQ. 1  ) 

M=1 

^  o  4  •  4  ■ 

IF 

(N.Ea.3) 

M=1 

555  •  4  • 

IF 

(N.EQ. 5) 

M  =  l 

556  a  ■  • 

IF 

(N.EQ. 2) 

M  =  0 

55  7  ... 

IF 

(N.EQ. 4) 

M  =  0 

5  53  ■  4  . 

IF 

(N.EQ. 6) 

M  =  0 

559  ■  .  . 

IF 

(N.EQ. 1  ) 

1=2 

5  6  0  ... 

IF 

(N.Ea.3) 

1  =  3 

.i  6  t  ... 

IT 

(N.EQ. 5) 

1  =  5 

5  C} ..  •  •  ■ 

IF 

(N.EQ. 2) 

1=5 

!  6  ?  ... 

IF 

(N.EQ. 4) 

1=5 

.76')  ... 

IF 

( N.EQ. 6) 

1  =  5 

5 1 5  ... 

IF 

'  N  .  E  Q  .  3  ) 

CALL 

.>6  6  ... 

IF 

(N.EQ. 5) 

CALL 

PLOTTING  OG  THE  SIX  GRAPHS 
CALL  Line  ( eei  ,aa. 1 3. i , m. i ) 
11  CONTINUE 


570 

571 
57^ 


CALL  PLOT  (6,5,0.0,999) 

RETURN 

END 
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